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TECHNICAL  NOTE  2604 

A  GENERAL  THEORY  OF  THREE-DIMENSIONAL  FLOW  IN 
SUBSONIC  AND  SUPERSONIC  TURBOMACHINES  OF 
AXIAL-,  RADIAL-,  AND  MIXED-FLOW  TYPES 
By  Chung -Hua  Wu 


SUMMARY 

A  general  theory  of  steady  three-dimensional  flow  of  a  nonviscous 
fluid  in  subsonic  and  supersonic  turbomachines  having  arbitrary  hub  and 
casing  shapes  and  a  finite  number  of  blades  is  presented.  The  solution 
of  the  three-dimensional  direct  and  inverse  problem  is  obtained  by 
investigating  an  appropriate  combination  of  flows  on  relative  stream 
surfaces  whose  intersections  with  a  z-plane  either  upstream  of  or  some¬ 
where  inside  the  blade  row  form  a  circular  arc  or  a  radial  line.  The 
equations  obtained  to  describe  the  fluid  flow  on  these  stream  surfaces 
show  clearly  the  several  approximations  involved  in  ordinary  two- 
dimensional  treatments.  They  also  lead  to  a  solution  of  the  three- 
dimensional  problem  in  a  mathematically  two-dimensional  manner  through 
iteration.  The  equation  of  continuity  is  combined  with  the  equation  of 
motion  in  either  the  tangential  or  the  radial  direction  through  the  use 
of  a  stream  function  defined  on  the  surface,  and  the  resulting  equation 
is  chosen  as  the  principal  equation  for  such  flows.  The  character  of 
this  equation  depends  on  the  relative  magnitude  of  the  local  velocity  of 
sound  and  a  certain  combination  of  velocity  components  of  the  fluid.  A 
general  method  to  solve  this  equation  by  both  hand  and  high-speed 
digital  machine  computations  when  the  equation  is  elliptic  or  hyperbolic 
is  described.  The  theory  is  applicable  to  both  irrotational  and  rota¬ 
tional  absolute  flow  at  the  inlet  of  the  blade  row  and  at  both  design 
and  off -design  operations. 


INTRODUCTION 

The  problem  of  three-dimensional  flow  in  turbomachines  of  axial-, 
radial-,  and  mixed-flow  types  is  treated  in  references  1  to  19.  Because 
of  the  enormous  mathematical  difficulties  involved  in  the  problem,  Lorenz 
(reference  l)  first  introduced  the  idea  of  an  infinite  number  of  blades 
of  infinitesimal  thickness  in  order  to  follow  the  flow  on  a  given  surface 
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Bauersfeld  (reference  2)  added  to  the  theory  the  condition  of  integra- 
bility  for  the  blade  surface  that  must  be  satisfied  in  the  inverse,  or 
design,  problem.  The  theory  is  further  clarified  and  strengthened  by 
the  works  of  Stodola  (reference  3),  von  Mises  (reference  4),  and  Dreyfus 
(reference  5) ,  and  is  the  basis  of  many  recent  investigations  on  axial-, 
radial-,  and  mixed-flow  compressors  and  turbines. 

For  incompressible  flow,  Ruden  (reference  6)  proves  that  the 
through-flow  solution  obtained  under  the  assumption  of  an  infinite  num¬ 
ber  of  blades  gives  a  circumferentially  average  value  of  the  fluid  prop¬ 
erties,  provided  the  deviations  of  the  fluid  properties  from  their  cir¬ 
cumferential  averages  are  small.  In  reference  7,  Traupel  points  out  the 
oscillatory  nature  of  radial  flow  in  a  multistage  turbomachine  and  gives 
solutions  of  the  three-dimensional  potential  flow  through  inclined 
stationary  blades  and  also  of  the  rotational  flow  through  a  homogeneous 
stage  of  identical  nontwisted  blades  for  an  incompressible  fluid  and  an 
infinite  number  of  blades  bounded  by  cylindircal  walls.  Meyer  gives  a 
detailed  treatment  of  three-dimensional  potential  flow  in  a  stationary 
blade  row,  for  an  incompressible  fluid  and  cylindrical  bounding  wall, 
in  reference  8,  where  the  solution  for  an  infinite  number  of  blades  is 
extended  to  a  finite  number  of  blades  by  the  vortex-and-source  method  of 
Ackeret,  which  is  originally  given  for  two-dimensional  flow  (reference  9). 
In  reference  10,  a  linearized  solution  for  an  incompressible  fluid  and  an 
infinite  number  of  blades  for  a  prescribed  loading  and  cylindrical 
bounding  walls  is  obtained  by  Marble,  and  is  used  later  to  investigate 
the  problem  of  mutual  interference  of  adjacent  blade  rows  and  off- 
design  operations  (reference  ll) .  Siestrunck  and  Fabri  (reference  12) 
also  obtained  a  linearized  solution  for  incompressible  flow,  and  the 
method  is  extended  to  compressible  flow.  For  general  wall  shapes, 
Spannhake  (reference  13)  examines  the  flow  through  diffuser  and  impeller 
by  the  use  of  bound  vortices  for  blades.  The  incompressible  through 
flow  in  a  mixed-flow  impeller  is  treated  by  Gravolos  (reference  14). 

In  reference  15,  Wislicenus  examines  the  influence  on  the  meridional 
flow  of  the  blade  force  and  nonuniform  circulation  along  the  blade  span. 

For  compressible  flow,  Reissner  (reference  16)  gives  a  blade-design 
method  in  which  the  extension  from  an  infinite  number  of  blades  to  a 
finite  number  of  blades  is  accomplished  by  the  use  of  a  power  series  in 
the  circumferential  direction,  and  the  terms  in  the  series  are  deter¬ 
mined  by  a  comparison  of  the  equations  for  an  infinite  number  of  blades 
and  a  finite  number  of  blades,  (in  reference  5,  Dreyfus  gives  a  method 
of  designing  water  turbines  of  thin  blades,  in  which  the  solution  for 
an  infinite  number  of  blades  is  extended  to  a  finite  number  of  blades  by 
the  use  of  a  power  series,  the  second  term  of  which  is  determined  from 
the  equations  of  continuity  and  irrotational  absolute  flow  and  is 
explicitly  given.)  In  reference  17  the  compressible  flow  problems  in 
axial  turbomachines  having  an  infinite  number  of  blades  are  treated,  and 
both  the  direct  and  inverse  problems’  are  considered.  Methods  for 
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limiting  solutions  for  zero  and  infinite  blade-row  aspect  ratios  and  a 
step-by-step  method  of  solution,  as  well  as  a  simpler  method  based  on  an 
approximate  knowledge  of  the  shape  of  the  streamline,  for  a  finite  blade - 
row  aspect  ratio  are  given.  Unaware  of  the  work  of  Traupel  at  the  time, 
the  authors  of  reference  17  also  emphasized  the  oscillatory  nature  of 
radial  flow  in  multistage  machines  and  suggested  the  use  of  a  simple 
sinusoidal  form  of  the  streamline  as  a  first  approximate  solution.  Their 
methods  are  derived  for  compressible  flow,  however,  and  are  also  extended 
to  the  case  where  both  the  hub  and  casing  walls  or  either  is  tapered. 
Reference  18.  gives  a  general  through-flOw  theory  for  both  direct  and 
inverse  problems  and  for  subsonic  or  supersonic  flow  in  turbomachines 
having  arbitrary  hub  and  casing  shapes.  The  supersonic  through  flow  in 
rotating  impellers  having  a  prescribed  flow  along  the  casing  and  pre¬ 
scribed  blade  shapes  is  treated  in  reference  19. 

A  general  theory  of  three-dimensional  flow  in  subsonic  and  super¬ 
sonic  turbomachines  of  axial-,  radial-;  and  mixed-flow  types  for  a 
finite  number  of  thick  blades  of  finite  thickness  has  been  developed  at 
the  MCA  Lewis  laboratory  and  is  presented  herein.  Both  the  direct  and 
inverse  problems  are  considered.  The  theory  is  applicable  to  either 
irrotational  or  rotational  absolute  flow  at  the  inlet  of  a  blade  row 
and  to  both  design  and  off-design  operations. 

In  the  section  BASIC  AEROTHERMODYNAMIC  RELATIONS,  the  motion  and 
energy  equations  for  the  unsteady  flow  of  a  nonviscous  compressible  fluid 
in  a  rotating  blade  row  are  expressed  in  terms  of  the  velocity  components 
and  of  two  basic  thermodynamic  properties  Of  the  fluid,  namely,  entropy 
and  a  modified  total  enthalpy  for  flow  in  rotating  blade  rows  with  change 
in  radial  distance  from  the  machine  axis.  Estimated  entropy  changes  due 
to  shock  waves  (in  the  case  of  supersonic  flow),  heat  transfer  (in  the 
case  of  a  cooled  turbine),  or  viscous  effect  can  be  easily  accomodated 
in  the  calculation.  The  equations  obtained  show  clearly  the  condition 
under  which  the  flow  through  blade  rows  can  be  treated  on  the  basis,  of 
irrotational  absolute  flow. 

In  the  following  section,  a  general  potential  equation  is  obtained 
for  steady  three-dimensional  compressible  flow  through  rotating  or 
stationary  blade  rows  when  the  absolute  flow  can  be  taken  as  irrota¬ 
tional.  The  methods  of  solution  for  both  subsonic  and  supersonic  flows 
are  briefly  discussed. 

A  simpler  method  of  solving  the  three-dimensional  irrotational 
(absolute)  flow,  which  is  also  applicable  to  rotational' (absolute)  flow, 
is  obtained  by  considering  fluid  flows  on  a  number  of  relative  stream 
surfaces  whose  intersection  with  a  z-plane  either  upstream  of  or  some¬ 
where  in  the  blade  row  form  a  circular  arc  or  a  radial  line.  Equations 
governing  the  flow  on  these  surfaces  are  obtained  in  the  next  four  sec¬ 
tions.  Through  the  use  of  a  stream  function  defined  on  the  stream 
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surface,  the  equations  of  continuity  and  motion  for  fluid  flow  on  these 
surfaces  are  combined  into  one  principal  equation.  The  character  of  the 
principal  equation  is  dependent  on  the  relative  magnitude  of  the  local 
velocity  of  sound  and  a  certain  combination  of  velocity  components. 

The  process  involved  in  solving  the  direct  and  inverse  problems  by 
this  approach  is  described  in  the  section  STEPS  FOR  COMPLETE  SOLUTIONS 
OF  THREE-DIMENSIONAL  DIRECT  AND  INVERSE  PROBLEMS.  In  the  inverse  prob¬ 
lem,  besides  the  blade  -thickness  distribution  determined  by  blade 
strength  and  other  considerations,  either  the  tangential  velocity,  a 
relation  between  the  tangential  and  axial  velocity,  or  one  other  rela¬ 
tion  is  prescribed  on  a  mean  stream  surface  about  midway  between  two 
blades.  The  last  section  gives  a  general  method  of  solution  of  the 
principal  equation  when  it  is  elliptic  or  hyperbolic. 


SYMBOLS 


The  following  symbols  are  used  in  this  report: 


a 


B,b 


rngi 

n  j 


C,c 


cp->  cv 

D_ 

Dt 

Dmq 


F,f 


H 

h 


velocity  of  sound 

integrating  factor  for  continuity  equation  for  S2  and  Sq 
surfaces,  respectively 

differentiation  coefficient  used  to  multiply  function  value 
at  point  j  to  give  the  m^  derivative  at  point  i  based 
on  n degree  polynomial 

nonzero  term  on  right-hand  side  of  continuity  equation  for 
S2  and  Sq  surfaces,  respectively 

specific  heat  of  gas  at  constant  pressure  and  volume, 
respectively 

differentiation  with  respect  to  time  following  relative 
motion  of  fluid  particle 

m^h  derivative  of  q 

vectors  having  the  unit  of  force  per  unit  mass  of  fluid 

given  function  of  ^u/^z  on  S2 

total  enthalpy  per  unit  mass  of  fluid,  h  +  ~V2 

static  enthalpy  per  unit  mass  of  fluid,  u  +  p/p 
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I  modified  total  enthalpy  for  flow  in  rotating  hlade  row  with 

change  in  radial  distance  from  machine  axis, 

h  +  i  o>2r2  or  H  -  oo(Vur) 

coefficients  of  first-  and  second-order  derivatives  in  the 
.  principal  equation 

k  thermal  conductivity 

L  distance  along  streamline 

l,  Cp  orthogonal  coordinates  on  surface  of  revolution 

M  mass  flow  between  mean  stream  surface  and  one  surface  of 

blade 

N  number  of  blades 

n  unit  vector  normal  to  relative  stream  surface  S 

p  static  pressure 

Q  heat  added  to  fluid  particle  along  its  path  per  unit  mass 

per  unit  time 

any  quantity  on  relative  stream  surface  S 
gas  constant 

remainder  term  of  m^1  derivative  at  point  i  obtained  by 
using  n^k  degree  polynomial 

r  radius  vector 

relative  stream  surface  passing  through  fluid  particles 

lying  on  a  circular  arc  upstream  of  or  midway  in  blade 
row 

S£  relative  stream  surface  passing  through  fluid  particles 

lying  on  a  radial  or  curved  line  upstream  of  or  midway  in 
blade  row 

entropy  per  unit  mass 
s/R 

static  temperature 


<1 

R 

“R1 

n 
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t  time 

U  velocity  vector  of  blade  element  at  radius  r 

u  interval  energy  per  unit  mass 

V  absolute  velocity  of  fluid 

W  velocity  of  fluid  relative  to  blade,  V  -  U 

V  Z^2  +  W^2 

x,y  independent  variables 

z  distance  along  turbomachine  axis 

.  1  Wu 

a  arc  tan  —  — * 

r  Wz 

X  ratio  of  specific  heats 

X  average  value  of  X  for  the  temperature  range  involved 

5  grid  spacing 

€  equal  to  1  and  r  for  and  S2  surfaces,  respectively 

£  independent  variable  z  or  r  for  surface  and  z  for 

Sg  surface 

T)  independent  variable  cp  and  r  for  S-j_  and  Sg  surfaces, 

respectively 

0  angular  distance  of  fluid  particle  measured  with  respect  to 

stationary  radial  line 

A  slope  of  characteristic  curves,  V  —3. 

A  tan  o 

.  a 

u  arc  sm  — 

w 

V  equal  to  r  and  1  for  and  Sg  surfaces,  respectively 

l  absolute  vorticity,  Vx  V 

p  fluid  density 
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^  generalized  variable  used  for  general  density  table 

cr  angle  between  tangent  of  streamline  or  boundary  wall  in  the 

meridional  plane  and  axial  direction 

T  radial,  axial,  or  angular  thickness  of  stream  sheet 

$  velocity  potential 

0  generalized  variable  used  for  general  density  table 

Cp  angular  distance  of  fluid  particle  measured  with  respect  to 
radial  line  on  rotating  blade 

X  angle  between  w  and  axial  direction 

stream  functions  defined  on  relative  stream  surfaces  S2  and  sl* 
respectively 

0)  angular  velocity  of  blade 
Subscripts: 
c  casing 

e  exit 

h  hub 

i  inlet 

l  meridional  component 

m  mean  stream  surface 

o  lower  limit  of  integration 

r,u,z  radial,  circumferential,  and  axial  components 
s  isentropic 

T  total  state 

T],£  components  in  tj-  and  ^-direction,  respectively 
1  on  S*l,  or  in  front  of  rotor 

on  S2, 


2 


or  behind  rotor 
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Superscripts: 

a,bj  .  .  .  k  refer  to  points  a,b,  .  .  .  k,  respectively 

BASIC  AEROTHERMODYNAMIC  RELATIONS1 

The  three-dimensional  flow  of  a  nonviscous,  compressible  fluid 
through  a  turbomachine  is  governed  by  the  following  set  of  basic  laws 
of  aerothermodynamics  •  From  the  principle  of  conservation  of  matter, 
the  equation  of  continuity  is 


|jr +V-(pW)  =  0 


(1) 


or 


V  •  W  +  -  jjg  P  =  0  (la) 

For  a  blade  rotating  at  a  constant  angular  velocity  co  about  the 
z-axis,  Newton* s  second  law  of  motion  gives 

^.co2r  +  2toXW  =  _  i Vp  (2) 

Because  the  boundary  walls  are  surfaces  of  revolution  and  the  relative 
flow  can  be  approximated  as  being  steady  in  many  cases,  it  is  convenient 
to  use  a  relative  cylindrical  coordinate  system  r,  cp,  and  z  with 
cp  measured  with  respect  to  the  rotating  blade  (see  fig,  1).  By  use  of 

5?  -  3?  +  (*-v>  *  '  3?  +  jvw2  -  wxCww) 


the  scalar  forms  of  the  equation  of  motion  (2)  in  the  axial,  circum¬ 
ferential,  and  radial  directions  can  be  expressed  as 

2 


r  awr  Wu  dWr  dWr  Wu"  ?  i  , 


W„  dw 


aw 


awu 

at” +  wr  a?- 


awu  wuawu  awu 
+  v  +  wz  ar 


wrwu  1 

+  2ctWr  =  - 
r  r 


_  dp 

pr  5V 


awz  awz 

at” +  wr  a?-  +  — 


aw, 


wu  ^z  TT  — z 

a^p” +  aa-  “ 


i  ap 

p  aa 


(2b) 

(2c) 


* 


^-Some  of  the  relations  given  in  this  section  have  been  given  in  refer¬ 
ence  18,  They  axe  repeated  here  for  completeness  and  easy  reference 
for  the  following  developments. 
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The  first  law  of  thermodynamics  may  be  written 


D(p ~1) 
Dt 


+  t)  _ L  =  o 

T\4-  ~  Jr  Th4-  ^ 


where  u  is  related  to  the  temperature  T  by 

du  ss  cv  dT  (4) 

and  Q  is  given  by  the  following  equation  if  only  conduction  is 
considered: 

Q  -  p-1  V-(kVT)  (5) 

For  the  ranges  of  temperature  and  pressure  encountered  in  ordinary 
turbomachines >  p,  p,  and  T  of  the  gas  are  accurately  related  by  the 
following  equation  of  state: 

p  =  R  pT  (6) 

Although  the  flow  of  the  gas  through  the  turbomachine  is  completely 
defined  by  the  preceding  equations  together  with  the  known  variations  of 
cv  and  k  with  temperature  and  the  given  boundary  and  initial  condi¬ 
tions,  it  is  found  more  convenient  in  references  17  and  18  to  express 
the  state  of  the  gas  in  terms  of  the  entropy  and  the  total  enthalpy  or 
a  quantity  I  of  the  gas,  besides  its  velocity  components.  These 
quantities  are  defined  as  follows:  . 


T  ds  a  du  +  p  d(p-1) 

H  =  h  +  |  V2 

I  =  h  +  i  W2  -  i  U2  =  H  -  o)(Vur) 


h  =  u  +  pp"1  (10) 

From  equations  (10),  (4),  and  (6)  is  obtained 

dh  =  (cv+R)  dT  =  cp  dT  =  dT  (11) 

where  y  is  equal  to  Cp/cy  and  is  a  function  of  temperature.  Another 
expression  for  dh  is  obtained  by  using  equations  (10)  and  (7),  so  that 


dh  =  ^  +  T  ds 


(11a) 
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By  the  use  of  equations  (7),  (4),  and  (6), 

1 


d(f)  =  7TTdlnP-7TidlnP  (12') 


and 


D  =  ^  d  m  T  -  d  m  p  (12a) 

can  he  obtained,  and  the  equation  of  continuity  can  be  -written 

V .  w  +  P  ■^n  y.  _  —  £  s-  o  ( 15) 

v  r-1  Dt  Dt  R  ' 

Equation  (13)  can  be  expressed  in  a  slightly  different  form.  From  the 
definition  of  the  local  velocity  of  sound  (reference  20 ) , 


a 


...  /Sp\ 

-  w 


By  the  use  of  equations  (12)  and  (6), 


a2  ts  x  ~ 
P 


(M) 


(14a) 


Substituting  this  relation  into  equations  (12a)  and  (13),  with  the  use  of 
equation  (ll),  results  in 


.  dh  s 

dlnp  =  ?-as 


(12b) 


and. 

.  W  +  _i_  _  JL  £  _  o 

V  W  a2  Dt  Dt  R 
From  equations  (9)  and.  (lla), 

-VP  +  JvW2  -  co2r  =  VI  -  Tv 

P  * 

The  equation  of  motion  (2)  can  then  be  written 

aw 


at 


-  WX(VXW)  +  2«XW  =  -  VI  +  Tvs 


(13a) 


(15) 


An  alternative  form  of  equation  ( 15) ,  which  involves  the  vorticity  of 
the  absolute  motion,  is  obtained  as  follows:  With  the  z-axis  parallel 
to  U) , 
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By  use  of  equations  (2),  (9),  and  (lla), 

DI  Ds  1  dp  1  ,  .  DW  DU 

Dt  =  TDt  +  pdt  +  p(WV)p+W’Dt-U*Dt 


“  T  §1  +  p  M  +W-(co2r-2a»XW)  -U-^f=Tg|  +  ±H||  (20) 


DU 


Ds  1  3] 


It  may  be  emphasized  that  the  preceding  equation  is  a  consequence 
of  the  equation  of  motion  (2)  and  the  thermodynamic  relations  (4),  (6), 
(7);  and  (10) .  For  steady  relative  flow,  the  rate  of  change  of  I 
along  the  streamline  is  seen  to  be  proportional  to  the  rate  of  change 
of  entropy  along  the  streamline. 

The  energy  equation  (3)  can  be  used  to  express  the  rate  of  change 
of  entropy  along  the  streamline  by  the  use  of  equation  (7)  as  follows: 

«  -  T  if  (21> 

The  preceding  equations  lead  to  several  important  general  consider¬ 
ations:  If  the  blade  rows  are  not  placed  too  close  together  and  no 
trailing  vortices  are  shed  from  preceding  blade' rows  (or  where  these 
effects  can  be  neglected) ,  the  fluid  properties  at  a  fixed  point  rela¬ 
tive  to  the  blade  can  be  taken  as  constant  with  respect  to  time.  Con¬ 
sequently,  according  to  equations  (20)  and  (2l),  the  quantities  s  and 
I  of  the  gas  remain  constant  along  the  streamline  for  adiabatic  flow. 
The  invariancy  of  I  means  that  the  rate  of  change  in  total  enthalpy 
along  the  streamline  is  equal  to  the  angular  speed  of  the  blade  multi¬ 
plied  by  the  rate  of  change  in  angular  momentum  (about  the  machine  axis) 
of  the  fluid  particle  along  its  streamline,  which  is  the  well-known 
Euler  turbine  equation  usually  derived  under  less  general  conditions. 

In  a  cooled  turbine  where  the  heat  transfer  may  be  large,  the  rate  of 
change  of  s  and  I  along  the  streamline  can  be  corrected  by  equa¬ 
tion  (2l)  for  an  estimated  value  of  Q.  Again,  for  steady  relative  flow, 
equation  (l5a)  shows  that  either  when  gradient  I  and  gradient  s  both 
vanish  or  when  the  difference  between  VI  and  TVs  vanishes,  the 
absolute  vorticity  either  vanishes  or  is  parallel  to  the  relative 
velocity. 

For  the  flow  through  a  stationary  blade  row  go  =  0,  W  becomes  V  , 

I  becomes  H,  and  equation  (l5a)  becomes 

-  VX  (V  xv)  a  -  VH  +  TV  s  (15b) 


which  agrees  with  similar  relations  previously  obtained  by  Vazsonyi 
(reference  21)  and  Hicks,  Guenther,  and  Vasserman  (reference  22).  It 
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is  interesting  to  see  that,  for  relative  flow  in  a  rotating  blade  row, 
VX(VXV)  becomes  WX(yXV)  and  H  becomes  I. 

If  it  is  assumed  that  the  fluid  enters  the  inlet  guide  vanes  of  a 
turbomachine  with  uniform  H  and  s  and  zero  vorticity  and  that  the 
flow  is  adiabatic,  s  does  not  vary  in  the  inlet  guide  Vanes  and  p 
is  then  a  function  of  only  p,  according  to  equation  (12).  Conse¬ 
quently,  by  virtue  of  Kelvin's  circulation  theorem,  the  absolute  vor¬ 
ticity  will  remain  zero  in  passing  through  the  inlet  guide  vanes  and  the 
flow  in  the  inlet  guide  vanes  can  be  treated  on  the  basis  of  irrota- 
tional  absolute  flow. 

If  the  guide  vanes  impart  a  radial  variation  of  tangential  velocity 
of  the  fluid  in  a  z -plane  downstream  of  the  vanes  similar  to  that  in  a 
potential  vortex,  that  is,  inversely  proportional  to  the  radius,  the 
circulation  is  constant  along  the  blade  span  and  the  fluid  maintains  a 
uniform  s  and  H  and  a  zero  vorticity  of  absolute  flow  entering  the 
following  rotor -blade  row.  If  the  rotor -blade  row  is  situated  far*  away 
from  the  inlet  guide  vanes,  the  fluid  enters  the  rotor  with  a  uniform 
I  'in  the  circumferential  direction,  as  well  as  in  the  radial  direction, 
and  the  flow  through  the  rotor  blades  can  again  be  treated  on  the  basis 
of  zero  absolute  vorticity  and  steady  relative  flow.  If  the  rotor  is 
close  to  the  guide  vanes,  however,  vortices  are  shed  from  the  inlet 
guide  vanes  because  of  periodic  variation  in  circulation  caused  by 
unsteady  flow,  and  the  flow  downstream  of  the  stator  and  through  the 
rotor  blades  should  theoretically  be  treated  on  the  basis  of  rotational 
flow. 

If  the  guide  vanes  impart  a  radial  variation  of  tangential  velocity 
of  the  fluid  at  a  z-plane  downstream  of  the  vanes  not  inversely  propor¬ 
tional  to  the  radius,  the  circulation  varies  along  the  span  of  the  guide 
vanes,  vortices  are  shed  from  the  trailing  edge  to  the  fluids  downstream 
in  the  direction  of  the  exit  velocity,  and  the  fluid  enters  the  follow¬ 
ing  rotor  blades  with  a  uniform  s  and  H  but  a  nonuniform  I  and  a 
nonzero  value  of  absolute  vorticity.  Consequently,  the  flow  through  the 
rotor -blade  row  can  no  longer  be  treated  on  the  basis  of  zero  absolute 
vorticity,  even  if  it  is  far  apart  from  the  preceding  guide  vanes. 

From  the  preceding  discussion,  the  choice  of  s  and  H  or  I 
as  the  two  basic  thermodynamic  variables  of  the  gas  besides  its  veloc¬ 
ity  components  is  apparent.  Compressor  and  turbine  rotors  are  usually 
designed  to  impart  or  subtract  the  same  amount  of  energy  to  or  from  the 
gas  radially;  hence  H  is  usually  radially  constant  throughout  the 
machine  if  the  inlet  flow  is  uniform  (except  in  the  boundary  layer  along 
hub  and  casing  walls).  If  the  circumferential  velocity  of  the  gas 
upstream  of  the  blade  row  is  zero  or  varies  inversely  with  radius,  I 
is  then  constant  throughout  the  machine.  These  facts  will  be  utilized 
in  the  following  developments. 
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POTENTIAL  EQUATION  FOR  THREE-DIMENSIONAL  FLOW  THROUGH 
ROTATING  BLADE  ROW 

Consider  first  the  special  case  of  steady  relative  flow  where  the 
fluid  upstream  of  the  “blade  row  is  free  of  vorticity  and  is  uniform  in 
I  and  s.  The  adiabatic  flow  through  the  blade  row  is  then  relatively 
steady  and  absolutely  irrotational  and  is  most  conveniently  treated  by 
the  use  of  a  velocity  potential  $  based  on  the  zero  absolute  vorticity 
and  related  to  the  relative  velocity  components  through  equation  (16)  as 
follows- 


SF 

1 


=  vr  =  wr 


r  5^"  Vu  "  wu  + 
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35  “  Vz  “  Wz 
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(22) 


For  steady  isentropic  flow,  the  continuity  equation  (l3a)  becomes 
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From  equations  (9)  and  (22), 

3®  1 


I  +  cp 


dcp 


/a®\2  ,  /l  3®\2  ^  /3®\2 

W_ 


3h 

35 


a2® 


i  a2®  wu  a2®  vu2  w 

c  r  ^r2  +  r  3r3<p  r  +  z  3555 


(23) 
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(25a) 
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By  the  use  of  equations  (22)  and  (25a)  to  (25c),  the  continuity  equa 
tion  (23)  may  be  written 
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(26) 


Equation  (26)  is  then  the  three-dimensional  potential  equation  for 
isentropic  .flow  in  a  rotating  blade  row0  It  is  seen  from  this  equation 
and  equation  (16)  that  the  three-dimensional  flow  through  a  rotating 
blade  row  cannot  be  treated  by  a  three-dimensional  flow  through  a  sim¬ 
ilar  stationary  blade  row  with  the  same  inlet  condition  relative  to  the 
blade  row,  as  in  the  case  of  two-dimensional  flow  on  a  cylindrical  sur¬ 
face,  because  the  difference  between  the  absolute  and  relative  vorticity 
2fc)  does  not  enter  into  the  two-dimensional  flow  on  a  cylindrical  surface 
but  does  enter  into  the  three-dimensional  case. 


Equation  (26)  is  very  similar  to  the  ordinary  three-dimensional 
potential  equation  for  flow  past  stationary  objects,  except  that  both 
relative  and  absolute  velocity  components  are  involved  in  the  coeffi¬ 
cients  of  $  derivatives  and  that  4>  is  directly  defined  by  the  abso¬ 
lute  velocity.  The  real  difficulty  in  solving  this  equation  lies  in 
the  fact  that  all  the  velocity  components  change  greatly  in  passing 
through  a  turbomachine  and,  consequently,  the  equation  cannot  be  lin¬ 
earized  and  yet  give  a  good  approximate  answer.  For  supersonic  relative 
velocity,  the  method  of  characteristic  surfaces  (references  23  and  24) 
may  be  used  to  solve  equation  (26),  with  the  initial  conditions  not 
given  on  a  characteristic  surface.  For  subsonic  relative  flow,  the 
equation  is  more  conveniently  written  in  the  form 


^  +  ±^+r<^+-ll+W  5  p 

dr2  +  r ■  5F  +  r2  ^2  +  Sz2  t  Wr 
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w  *  .. 

u  o  In  p 
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+  wz 


h|-p » o 


(26a) 


and  can  be  solved  by  Southwell's  relaxation  method  (reference  25)  or 
other  numerical  methods  using  the  differentiation  formulas  obtained  in 
reference  26  to  take  care  of  the  unequal  grid  spacings  near  the  blade 
surfaces  and  the  curved  hub  and  casing  walls.  The  last  three  terms  in 
equation  (26)  are  computed  from  the  $  values  or  velocities  obtained 
in  the  previous  cycle  and  kept  as  constants  during  each  improvement  of 
4>  values,  and  the  whole  process  is  repeated  until  the  desired  accuracy 
is  obtained.  Because  a  three-dimensional  stream  function  cannot  be 
defined,  the  use  of  velocity  potential  results  in  a  boundary- value  prob¬ 
lem  of  the  second  kind,  which  is  more  difficult  to  handle  in  the  calcu¬ 
lation  than  the  first  kind.  The  boundary  condition  to  be  satisfied  is 
that  the  relative  velocity  normal  to  the  moving  blade  is  zero,  or 
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Wr  nr  +  Wu  nu  +*WZ  nz  *  0  (27) 

where  n  is  the  unit  normal  vector  at  the  boundary  surface ,  and  that, 
at  inlet  and  exit  stations  far  away  from  the  blade,  the  velocity  is 
parallel  to  the  bounding  hub  and  casing  walls,  which,  in  the  case  of 
the  axial  machine,  means  that 

Vp  =  ||  =  0  at  z  =  i  oo  (27a) 

In  both  the  subsonic  and  supersonic  cases,  the  solution  is 
extremely  time-consuming.  Furthermore,  this  direct  approach  to  the 
three-dimensional  problem  requires  that  the  absolute  velocity  at  the 
inlet  to  the  blade  row  be  irrotational  and  the  flow  be  adiabatic.  In 
actual  machines,  the  flow  entering  the  blade  rows  is  always  rota¬ 
tional,  which  is  caused  by  a  nonuniform  total  enthalpy  and  entropy  at 
the  inlet  of  the  machine,  by  entropy  change  caused  by  shock  waves  or 
heat  transfer,  or  by  the  effect  of  boundary  layers  along  the  hub  and 
casing  walls.  Some  other  approach  to  the  problem,  which  is  simpler  to 
handle  and  is  also  applicable  to  rotational  inlet  flow,  is  therefore 
desirable.  One  approach  is  suggested  in  the  following  sections. 


FOLLOWING  FLUID  FLOW  ON  RELATIVE  STREAM  SURFACES 

In  order  to  solve  the  steady  three-dimensional  flow,  with  either  an 
irrotational  or  rotational  absolute  flow  at  the  inlet,  in  a  relatively 
simple  manner,  an  approach  is  taken  to  obtain  the  three-dimensional 
solution  by  an  appropriate  combination  of  mathematically  two-dimensional 
flows  on  essentially  two  different  kinds  of  relative  stream  surface 
(figs.  1  to  3) .  The  first  kind  of  relative  stream  surface  is  one  whose 
intersection  with  a  z-plane'  either  upstream  of  the  blade  row  or  midway 
in  the  blade  row  forms  a  circular  arc  (fig.  l)  .  The  second  kind  of 
relative  stream  surface  is  one  whose  intersection  with  a  z-plane  either 
upstream  of  the  blade  row  or  somewhere  inside  the  blade  row  forms  a 
radial  line  (fig.  2).  These  two  kinds  of  relative  stream  surface  will 
be  hereinafter  designated  stream  surfaces  Sp  and  S2j  respectively. 


Sj_  Stream  Surface  of  First  Kind 

In  figure  1  is  shown  a  stream  surface  of  the  first  kind  formed  by 
fluid  particles  lying  on  a  circular  arc  ab  of  radius  oa  upstream 
of  the  blade  row.  It  is  usually  assumed  in  ordinary  two-dimensional 
treatments  (for  example,  references  27  to  30)  that  the  stream  surface 
thus  formed  is  a  surface  of  revolution.  In  the  following  development, 
the  surface  will  be  allowed  to  take  whatever  shape  it  should  have  in 
order  to  satisfy  all  the  equations  governing  the  three-dimensional  flow. 
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In  most  cases ,  the  deviation  of  the  surface  from  a  surface  of  revolution 
is  not  large,  and  it  is  satisfactory  to  consider  Sq  surfaces  formed  by 
fluid  particles  originally  lying  on  a  circular  arc  upstream  of  the  blade 
row.  If  the  rotationality  of  the  inlet  absolute  flow  is  large,  if  the 
blade  is  designed  for  a  velocity  diagram  quite  different  from  the  free- 
vortex  type,  or  if  the  blade  length  is  long  in  the  direction  of  the 
through  flow  (radial-  and  mixed-flow  machines),  the  twist  of  the  surface 
maybe  quite  large,  resulting  in  very  large  circumferential  derivatives . 
If  this  effect  is  found  during  calculation  or  known  from  experience,  it 
is  more  satisfactory  to  consider  Sq  surfaces  formed  by  fluid  particles 
originally  lying,  in  front  of  the  blade  row,  on  curves  inclined  to  the 
circular  arc  in  a  direction  opposite  to  the  twist  of  the  surface.  In 
this  way,  the  intersection  of  the  Sq  surface  with  a  constant  z-plane 
about  midway  in  the  flow  path  is  a 'nearly  circular  arc,  and  the  total 
twist  of  the  surface  will  be  about  equally  distributed  toward  the 
upstream  and  downstream  directions  (fig.  3).  If  this  distribution  of 
the  twist  of  the  stream  surface  is  still  not  enough,  it  may  be  necessary 
to  divide  the  complete  flow  path  into  a  few  shorter  paths  and  consider 
an  Sq  surface  for  each  of  them.  Under  these  conditions,  Sq  surfaces 
formed  by  fluid  particles  originally  lying  on  the  hub  or  casing  walls 
upstream  of  the  blade  row  should  not  be  chosen  in  order  that  the  compli¬ 
cation  arising  from  the  possibility  of  fluid  particles  leaving  the  wall 
and  flowing  along  the  blade  surface  may  be  avoided.  In  such  cases  it 
is  better  to  consider  the  Sq  surface  a  short  distance  from  the  hub  and 
casing;  otherwise,  for  an  approximate  solution  the  fluid  can  be  con¬ 
sidered  to  follow  the  hub  and  casing  walls,  which  are  surfaces  of  rev¬ 
olution,  and  the  calculation  is  thus  much  simpler  than  that  for  a 
general  surface. 


S2  Stream  Surface  of  Second  Kind 

A  stream  surface  of  the  second  kind  is  shown  in  figure  2.  The 
most  important  surface  of  this  family  is  the  one  about  midway  between 
two  blades  dividing  the  mass  flow  in  the  channel  into  two  approximately 
equal  parts.  This  surface  is  designated  the  mean  stream  surface 
(S2,m)  •  For  blades  with  radial  elements,  such  as  the  one  shown  in  fig¬ 
ure  2,  it  is  convenient  to  consider  a  mean  stream  surface  formed  by 
fluid  particles  originally  lying  on  a  radial  line  ab  upstream  of  the 
blade  row  if  the  twist  of  the  surface  is  not  expected  to  be  large. 
Otherwise,  the  radial  line  is  chosen  about  midway  in  the  passage  with 
the  fluid  particles  originally  starting  out  from  a  curved  line  upstream 
of.  the  blade  row  such  as  shown  in  figure  3. 

The  mean  stream  surfaces  for  axial-flow  gas  turbines  designed  on  a 
free -vortex  velocity  diagram  are  shown  in  figures  3  and  4.  The  radial 
element  of  the  mean  stream  surface  (fig.  4)  is  chosen  accordingly  as  the 
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stator  is  designed  to  aline  the  blade  sections  radially  at  the  leading 
edge,  trailing  edge,  or  somewhere  between.  Inasmuch  as  the  rotor -blade 
sections  are  usually  alined  radially  at  or  near  the  center  of  gravity 
of  the  blade  sections,  the  radial  position  of  the  mean  relative  stream 
surface  is  chosen  at  the  same  position  (figs.  3  and  4).  The  continua¬ 
tion  of  the  stream  surface  outside  the  blade  row  is  not  shown.  The  mean 
stream  surfaces  for  the  inlet  stage  of  a  multistage  axial  compressor 
designed  on  the  principle  of  a  symmetrical  velocity  diagram  at  all  radii 
are  shown  in  figure  5. 

Both  of  these  two  kinds  of  stream  surface  are  employed,  in  general, 
in  the  solution  of  the  three-dimensional  problem.  The  correct  solution 
of  one  surface  often  requires  some  data  obtainable  from  the  other,  and, 
consequently,  successive  solutions  between  these  two  are  involved.  Yet, 
the  solution  of  each  surface  is  manageable  with  the  present  mathematical 
technique  and  computational  facilities.  In  many  practical  cases,  and 
especially  in  the  inverse  problem,  however,  this  iteration  may  not  be 
required  if  only  an  approximate  solution  is  required  or  if  the  prescribed 
values  lead  to  a  satisfactory  blade  shape.  These  points  will  be  dis¬ 
cussed  in  the  section  next  to  the  last  (pp.  53  to  57), 


Relations  among  Relative  Velocity  of  Fluid,  Coordinates  of 

Stream  Surface,  and  Normal  to  Stream  Surface 

In  general,  the  coordinates  of  the  stream  surfaces  and  , their  differ¬ 
entials  are  related,  respectively,  by  the  following  equations: 

S(r,cp, z)  a  0  (28) 

^dx  +  Hdqp+Mdz“°  ^ 

Rather  than  use  the  three  partial  derivatives  of  S  with  respect  to 
the  coordinates,  it  is  convenient  to  consider  the  unit  vector  n  normal 
to  the  surface,  which  is  related  to  S  by 

“r  _  nu  _  nz 

3s  “  I“3s  “  ciS  “  r — o - — 

s  s  a|(!)  ♦(£ 

The  vector  n  is,  of  course,  perpendicular  to  the  relative  velocity  W, 
so  that 


n  .W  *  0 
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or 


nr^r  +  nu^u  nz^z  “  0  (3l) 

By  using  equation  (30),  equation  (29)  can  "be  written 

Up  dr  +  nu  r  dcp  +  nz  dz  =  0  (29a) 

The  vectors  n  and  W  are  shown  on  and  S2  surfaces  in  fig¬ 

ures  1  and  2. 


EQUATIONS  GOVERNING  FLUID  FLOW  ON  8±  SURFACE  FOR  AXIAL-FLOW  AND 
AXIAL-DISCHARGE  MIXED-FLOW  TURBOMACHINES 

If  the  fluid  motion  on  8-^  is  followed,  equations  (28)  and  (31) 
can  be  used  to  eliminate  one  of  the  three  coordinates.  For  axial-flow 
(figs.  6(a)  to  6(c))  and  axial-discharge  mixed-flow  turbomachines 
(fig.  6(d)),  it  is  convenient  to  express  r  in  terms  of  Cp  and  z. 

For  radial-flow  and  radial-discharge  mixed-flow  turbomachines  (figs.  6(e) 
and  6(f)),  this  system  will  encounter  difficulty  at  the  exit  where  the 
rate  of  change  of  fluid  state  with  respect  to  z  becomes  infinite.  It 
is  therefore  necessary  to  eliminate  z  and  to  consider  r  and  cp  as 
the  two  independent  variables. 


Flow  Along  General  8±  Surface 

For  axial-  and  mixed-flow  turbomachines,  any  quantity  q  on  the 
S]_  surface  is  considered  a  function  of  cp  and  z;  that  is, 

q.  =  <l[<ftz,r(cp,z)[] 

The  change  in  q  along  due  to  a  small  change  in  op  while  z  is 

held  constant  is  (see  fig.  l) 

d<1  =  ||  acp  +  ^  ||  aqp 

From  equation  (29a), 
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hence 


dq. 

Similarly,  for  dcp  =  0, 

dq 


/Sq_  V  Sq\ 

\Scp  ~  rij,  Sr / 


dcp 


/Sq  nz  Sq 

\Sz  nr  Sr 


dz 


With  a  bold  partial  derivative  sign  used  to  denote  the  rate  of 
change  of  any  quantity  q  on  S]_  with  respect  to  cp  or  z ,  with  the 
other  kept  constant,  the  preceding  relations  give 


1  3q  1  dq  nu  Sq 

r  acp=  r  S<p  ~  n^.  Sr 

3q  _  Sq  nz  Sq 
3  z  ~  ciz  Sr 


~\ 


J 


(32) 


With  the  relations  (3l)  and  (32),  the  rate  of  change  of  q  along  a 
streamline  on  S-^  is 


2a  _  4.  w  2a 

Dt  "  r  acp  z  az 


(33) 


Equations  of  continuity  and  motion.  -  When  the  fluid  motion  is 
followed  along  the  stream  surface  and  equations  (31)  and  (32)  are  used, 
the  continuity  equation  for  steady  relative  motion  becomes 


where 


c(«P. 


a(Pwu) 

3(PWZ) 

acp  + 

8  z 

S(Wrr) 

\r 

dr 

=  P  c(cp,z) 


Swu  Sw\ 

+  nu  3F”  +  nz  3F7 


(34) 


(35) 


For  rotational  steady  relative  motion,  the  equations  of  motion  (14) 
in  the  radial,  circumferential,  and  axial  directions  are 
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<  wuSw  /S wr  Swz\  w  Si  _ 

•  wu  ST  +  T  w  z  W"  “  3F7  "  u  =  "  5?  +  T 


\5S~  “  3F~, 


Swr  Swz' 


jS  +  l7 

or  or 


wrwu 


oWu  wr  Swr 


r  Sr  r  <^cp 


1  Swz  Sw, 


i  Si  t  Ss 


w  “  Try  +  a^r  -  -  7  ^  +  F 


&.^+„  + 
r  \6z  Sr  /  u  \r  S<P  dz  /  3z 


Relations  (9),  (16),  (3l),  and  (32)  along  the  relative  stream  sur¬ 
face  S]_  can  he  used  to  reduce  equations  (36)  to  the  following: 


Wu  Wu  3Wp 

-  .j.  —  —  +  W7  ~z - 2a)Wn  = 

r  r  3cp  z  3z  u 


Sh  Ss  2 
Sr  -  T  ST”  r. 


wrwu  wr  wr 

r  ’  r  9cp 


z  \r  3Cp  "  3z  , 


TT  I  -  -  -I  .  <5  TT  1  31  T  3s  nu  /Sh  Ss  2,, 

_Wz  >  +  2wWr  ='F  3^+?a^“5;  VSF"t  SF"“r 


wr  37-  +  w, 


l^z  fu’ 

r  3<P  3z  , 


z  /Sh  _  SS  ...2, 


«  -£=  +  T^-  — 

3z  8z  \or  or  y 


The  last  term  in  each  of  the  preceding  three  equations  is  propor¬ 
tional  to  the  components  of  the  normal  vector  and  therefore  can  he 
expressed  as  a  component  of  a  vector  that  is  parallel  to  n  and  has 
the  dimension  of  force  per  unit  mass.  If  this  term  is  defined  as 


1  /Sh  _  Ss  2 


r)  n  o  -  -L/i 
/  %  \P 


the  preceding  equations  can  he  written 


wuaw  awr  w* 

~  3Cp  +  Wz  57"  "  r  ^u  "  fr 


(39a) 


WrWu  ¥r  3Wr  ^  /  x  3WZ  3WU\  _  1  31  .  T  3s  .. 


r  3cp 


r  9cp  3z 


+  2oiWr  *  -  ±  S  +  i  |r*+  f„  (39b) 


37  +  wu  {■ 


if^z  _  !u 
r  3cp  az  . 


r  3cp  r  sep  u 


|i  +  T  |£  +  f 

3z  3z  z 


(39c) 


22 


'MCA  TN  2604 


Similarly,  the  equations  of  motion  in  the  form  of  equations  (2)  can  he 
written 


where 


?k  Ll_,,  > 

r  acp  z  Qz  r  “  r  ] 


Wu  3Wu 
r  acp 


+  Wz 


3WU  WrWu 
az  +  ~r 


2aSWr 


l  ap 

pr  acp 


+  fu 


wu  awz  awz 
7  5T  +  wz3T 


i  ap 
p  az 


+  f 1 
T  iz 


J 


(40) 


ii|En 

%  p  or 


Because  this  vector  f  is  parallel  to  n ,  it  is  perpendicular  to 
the  relative  velocity  of  fluid,  or 


frWr  +  fuWu  +  fzWz  -  0 


(41) 


By  the  use  of  equation  (41)  and  equations  (39),  it  can  he  shown  that  for 
steady  flow  on  an  S-]_  surface. 


/ 


M  _  Ds 
Dt  ~  1  Dt 


(41a) 


which  agrees  with  equation  (20)  .  Therefore,  for  the  present  problem  of 
steady  relative  flow  on  a  stream  surface,  the  relation  (41a)  can  be 
taken  either  as  one  of  the  equations  of  motion  or  to  represent  the  rela¬ 
tion  given  by  equation  (41)  .  In  other  words,  there  are  only  four  inde¬ 
pendent  relations  among  equations  (39a),  (39b),  (39c),  (41),  and  (41a). 

Just  as  in  the  case  of  the  continuity  equation,  either  set  of  the 
preceding  equations  of  motion  is  expressed  in  terms  of  the  special  par¬ 
tial  derivatives  with  respect  to  the  two  independent  variables  Cp  and 
z.  The  effect  of  radial  pressure  gradient  is  taken  into  account  in  all 
these  equations  by  the  f  term,  which  is  neglected  in  the  ordinary  two- 
dimensional  treatment  on  a  surface  of  revolution.  Equations  (28),  (31), 
(34) ,  and  (39)  or  (40) ,  however,  lead  to  a  possibility  of  correctly 
solving  the  three-dimensional  flow  of  fluid  particles  on  an  S]_  surface 
in  a  mathematically  two-dimensional  manner. 


Principal  equation.  -  The  equations  of  continuity  and  the  equation 
of  motion  in  the  circumferential  direction  can  be  combined  into  a  prin¬ 
cipal  equation  through  the  use  of  a  stream  function  \|r  as  follows: 
First,  if  a  variable  b  is  introduced  such  that 


A 


\ 
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D  In  b 


nz  Wz 
C  +  ~r" 


nz  M 
C  "  “r  r/ 


if  nz  wz\ 

W  V  "  Dr  r  / 


dx  (42a) 


in  which  the  integration  is  performed  along  a  streamline  on  the  Sq_ 
surface,  then  the  continuity  equation  (34)  can  be  written 


3(bpW,  )  3(bpW  r) 

- H-  + - o 

acp  az 


(34a) 


The  preceding  equation  is  the  necessary  and  sufficient  condition  that 
there  exist  a  function  \|r  with 


§5«  ^z 


(43a) 


H  -  -bP»u 


(43b) 


The  difference  in 


at  two  points 


on  the 


surface  is 


hp(Wzr  hep 


In  particular,  the  difference  in  i|r  at  two  points  j  and  k  on  the 
constant-z  plane  at  the  inlet  where  the  fluid  state  is  uniform  is 


-  t'5  -  *>iPi  WZ;i  r  dcp 

Jcpj 

These  two  equations  show  that,  physically,  the  integrating  factor  h 
can  be  interpreted  as  proportional  to  the  local  radial  thickness  of  a 
thin  stream  sheet  whose  mean  surface  is  the  stream  surface  considered 
here.  The  continuity  equation  (34a)  can  also  be  obtained  by  consider¬ 
ing  the  mass  flow  going  into  an  element  of  such  a  stream  sheet  as  shown 
in  figure  7.  By  equating  to  zero  the  mass  flow  going  into  the  element, 
which  is  defined  by  two  axial  planes  dcp  apart  and  two  normal  planes 
dz  apart  (see  fig., 7(a)),  and  letting  dcp  and  dz  approach  zero, 
there  is  obtained  , 
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8(tpWu)  3(t  pW  zr ) 

9Cp  +  9z 


(34b) 


■where  T  is  the  radial  thickness  of  the  stream  sheet.  From  equa¬ 
tions  (34a)  and  (34b),  it  is  apparent  that  b  is  proportional  to  T, 
and  the  differences  in  ijr  at  two  points  j  and  k  as  given  by  the 
two  equations  preceding  equation  (34b)  are  proportional  to  the  mass  flow 
across  any  line  joining  the  two  points,  Jn  actual  computation,  only  the 
ratio  b  to  b^  or  T  to  is  important  (a  different  initial  value 

amounts  to  a  different  constant  multiplier  of  the  relation  between  \|r 
and  mass  flow)  •  In  the  following,  b  will  be  retained  in  the  equation, 
but  in  actual  calculation  it  is  simpler  to  evaluate  the  ratio  T  to 
than  to  evaluate  the  ratio  b  to  b^,  both  from  the  data  obtained 

on  the  Sr,  surface  to  be  discussed  later.  Although  the  evaluation  of 

this  ratio  requires,  in  general,  calculations  on  the  S2  surfaces,  a 
means  is  nevertheless  provided  to  determine  correctly  the  flow  on  a 
general  S-^  surface  through  iteration. 

From  equation  (43), 


•h  1  a¥z  _  18 1  8  In  bp  3\|f  nu  1  3\|r 

r  8Cp  r2  ac(£  '  r2  8Cp  9cp+  r2  8cp 

,  9¥u  _  3^  3  In  bp  3\|r 

^  3z  ”  _  2  3z  '  8z 

9z 


(44a) 


(44b) 


The  third  terms  in  the  preceding  equations  can  be  expressed  in  terms 
of  h  through  the  use  of  equation  (l2b): 

dlnbp=dlnb+idh-ds*  (45) 

a  c 

where  s*  =  s/R.  But  from  equations  (9)  and  (43), 


h  a  I  + 


m2r2  Wr‘ 


2  2 

Then  from  equations  (45)  and  (46), 


-!(M-2[(if!)2  +  (!i)2]  (-> 


[*2  -  (wu 


2  +  wz2 


Jj  1  3  In  bp  _  1  ^I+ 


3Cp 


r  3cp 


o>2r2  "  Wr2\  a2 
- 2 - 1  +~ 


W  2  n, 


/a  In  b  3s*\ 

\  3Cp  ‘  3cp  j 


/  v-l[Wz32t  Wu92f  )  "z  iXu 

^  ^  Vr2  3<?£  "  r  "  r  “r 


U 
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Substituting  the  preceding  two  equations  into  equations  (44a)  and  (44b) 
and  adding  yield 


(a2-W  2)  24- 


2 

+  w 

2\~]  /  1  3Wz 

Wu\ 

a2  -  V  a2^ 

,  WuWz  a2 

u 

z  /J  \r  9Cp 

9z  / 

2 

r 

acp2 

r  az 

1 

9 

I  m2r2-W  2N 

II  x  i  £*  . 

Ib  In  b 

8s  *\ 

^-Wu2  % 

r 

8cp 

\I+  2  / 

\  9cp 

8Cp/  _ 

r  nj. 

cn2r2  - 


8z 


X  + 


2  /  8  In  b  3s*\  ^z  nz 
/  +  a  \  3Z  -9z/)~.r  n^ 


r  acp 


-H  <«> 


Substituting  the  preceding  equation  into  equation  (39b)  and  dividing  by 
a2  give  the  principal  equation  for  the  determination  of  fluid  motion 
along  a  general  surface: 


1  - 


wuwz  a2*  (  Wz  )  a2*  N  at  M  8Jr  „  n 


■where 


(48) 


w  9  In  b  ,  9s 
^  "  3z  9z 


*  .  1  (  31  ,  v 

"  +  3  V"  +  Wr  9V 


cjb2r2  +  Wz2  nz 


1  a  In  b  ,  1  8S*  ,  1  f  1  91  9Wr  a2-(Wu2+Wz2)  Wr 
W="r  acp  +  r  9Cp  +a2  \  r  acp  r  9cp  r  Wz 


a 


_(WU2-Kb2r2)nu\  a 


5 zr 


\  2 

■  wz 


r_  1  91  .  T  9s  -  /  l  9Wr 

|_  r3Cp  r  8cp  u+r^r  3cp 


The  equation  of  the  characteristics  of  the  differential  equation  (48)  is 
(reference  31) 
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Equation  (50)  shows  that  the  characteristics  are  real  when 
a/wu2  +  wz2  >  a,  in  which  case  the  method  of  characteristics  for  two  inde¬ 
pendent  variables  (references  20,  30,  31^,  and  32)  can  be  applied.  When 


/\|wu2  +  WZ2<  a. 


the  characteristics  are  imaginary,  and  it  is  more  con¬ 


venient  to  solve  the  equation  by  relaxation  (references  25,  33,  26,  and 
29)  and  matrix  methods  (references  26  and  29)  in  the  following  form, 
which  is  obtained  by  substituting  equations  (44)  into  equation  (39b): 


1  32i|f  a2^  [~1  /a  In  bp  nu\  3\|r  3  In  bp  3^1 

r2  8cp2  3z2  "  Lr2  \  acP  "  V  acp  +  dz  3zJ  + 

(bp)2 
i  at 
r  acp 


1 

r 


H  +  **L  +  f  + 

acp  +  r  acp  +  u  + 


wr  awr  /wu  \ 
7a|-wri,7  +  2ffijJa0  (48a) 


Procedure  of  solution,  -  It  may  be  noted  that  equation  (39b), 
instead  of  (39c),  is  chosen  to  form  the  principal  equation  (48)  or  (48a), 
because  fu  is,  in  general,  much  smaller  than  fz.  The  various  quan¬ 
tities  appearing  in  equation  (48)  or  (48a)  are  to  be  computed  from 
other  equations  given  earlier .  With  the  introduction  of  the  stream 
function,  there  are  altogether  seven  basic  independent  relations  -  one 
energy  equation  (21) ;  three  equations  of  motion,  (48)  or  (48a),  (39a), 
and  (39c) ;  two  equations  between  \|r  derivatives  and  fluid  properties, 
(43a)  and  (43b);  and  the  orthogonal  relation  between  W  and  f,  equa¬ 
tion  (41)  or  (41a) .  On  the  other  hand,  there  are  ten  basic  dependent 
variables  in  \| r,  b,  Wr,  Wu,  Wz,  fr,  fu,  fz,  s,  and  I  (or  h  '  or  p) 

to  define  the  flow  and  the  shape  of  the  surface.  In  general,  the  vari¬ 
able  b  is  to  be  evaluated  according  to  equation  (42a)  or  from  the  var¬ 
iation  in  the  radial  thickness  of  stream  sheet  using  the  data  obtained 
in  the  solution  of  S2  surfaces  and  is  therefore  considered  as  given 
here.  If  during  the  complete  solution  of  the  three-dimensional  flow  the 
shape  of  an  S]_  surface  is  taken  as  the  one  obtained  by  joining  corre¬ 
sponding  streamlines  obtained  on  Sg  surfaces  of  the  preceding  cycle, 
two  relations  between  the  n-  (or  f-)  components  are  given  by  equa¬ 
tions  (29),  and  there  are  now  altogether  nine  equations  to  be  solved  to 
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find  the  nine  unknowns.  Alternatively,  the  variation  of  Wr  may  be 
considered  as  known  from  the  Sg  solutions  of  the  previous  cycle,  and 
the  remaining  eight  variables,  which  determine  the  flow  on  and  the  shape 
of  the  S-|_  surfaces,  can  be  determined  from  the  seven  preceding  rela¬ 
tions  given  and  the  following  additional  relation:  Because  fr,  fu, 

and  fz,  respectively,  are  proportional  to  i  and  ^  of  the 

integral  surface  S,  they  satisfy  the  following  equation  (reference  34) 

f.VXfaO  (51) 

which  may  be  written 


1  ^z 

1  S(fur) 

5fz\  _ 

1  S(fur) 

_  i  ^fr 

_r  Sep 

r  Sr  _ 

+  fu  - 

Sr  /  +  fz 

jc  dr 

~  r  S<P_ 

By  using  equations  (31)  and  (32),  equation  (51a)  becomes  simply 

a  JL  (fur'\ 

3cp\f r/  3z  \fr  / 


(51a) 


(51b) 


This  equation  can  he  used  to  give  fu  hy  integrating  along  a  constant 
cp  line : 


If  at  z  -  zQ,  f u  «  0,  then 

nz 


(51d) 


In  this  case,  then,  the  shape  of  the  surface  is  determined  after 

the  f -components  (or  n- components)  are  obtained  in  the  solution.  In 
either  case,  equations  (21)  and  (41a)  are  invariably  to  be  used  first  to 
determine  the  change  of  s  and  I.  If  the  flow  is  isentropic,  s 
and  I  remain  constant  along  its  streamlines  on  the  surface.  (For 
such  a  case  and  for  a  uniform  inlet  condition,  p  in  the  continuity 

equation  may  be  replaced  by  h'  and,  consequently,  the  <P  and  z 
derivatives  of  s,  as  well  as  those  of  I,  will  not  be  involved  in 
the  equations  (46)  to  (48) ) .  In  case  of  heat  transfer  or  shock,  the 
changes  in  s  and  I  can  be  estimated  by  whatever  method  is  available 
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and  used  in  the  calculation.  For  supersonic  flow,  all  the  equations  are 
used  to  compute  the  fluid  state  at  each  point  and  the  solution  is  carried 
downstream  step  by  step.  For  subsonic  flow,  iteration  over  the  whole 
domain  is  necessary.  The  details  of  these  computations  will  he  given  in 
the  last  section.  In  general,  the  solution  of  the  flow  on  the  general 
Sj_  surface  is  very  laborious,  and  is  to  be  used  in  the  final  stages  of 
calculation  of  the  complete  three-dimensional  problem  or  when  a  high¬ 
speed  computing  machine  is  available. 

If  the  flow  is  such  that  it  may  be  assumed  to  take  place  on  a  sur¬ 
face  of  revolution  (at  the  hub  and  casing  walls  or  other  radii),  the 
equations  are  considerably  simplified  as  follows: 


Flow  along  Surface  of  Revolution 


When  the  S]_  surface  is  a  surface  of  revolution. 


n 


u 


-  fu  = 


0 


(52) 


Let 


nz 

Hr 


tan  o 


A 


(53) 


where  X  is  a  given  function  of  z.  (For  a  conical  flow  surface,  X 
is  simply  a  constant.)  Equation  (35)  now  gives 


Whether  c  can  be  taken  as  zero  will  be  determined  by  the  relative 
magnitude  of  the  three  terms  on  the  right  side  of  the  equation.  In 
general,  for  nonnegligible  c,  equations  (43)  now  become 


"5cp 


rbpWz 


(55a) 


at 

at  -  -  bew 


(55b) 


Because  Wr  is  now  related  to  Wz  by  equation  (53),  the  three  velocity 

components  can  be  solved  simultaneously  as  follows:  By  use  of  the  rela¬ 
tions  (52)  and  (53),  equation  (39b)  can  be  changed  to 


(l+X2) 


1 

r  5qT 


8W 

dZ 


u  \  /wu  L  „  \  1  /l  dl  T  ds\  .  ,ccs 

-  -  ^  V7  3?  '  ?  3W  "  0  (S6) 


990?. 
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Instead  of  equation  (46) , 


i  2 

h  =  I  +  — - tt  (tp) 


should  he  written.  Then 


1  Si  ,  a 2  /dint  Ss*\  .  s-lfl  ,\2^  Wz  ^2\|r  Ww  a^t_ 

7  3cp+-  \r^~'WJ  "  (bp)  (1+  }  3  ^ 


,  2  T,2X  1  S  In  bp  1  5l  ,  a^  h  In  ti  Ss 
(a  -W  )  -  -  r  r  y  acp  ~  3q 


(a2-«2>  *£*-§♦**♦.*  (±£*  -|f)  +  («/«z2)  ^WA  §  - 


\  -1  /  ,  \2  \  *z  9  \|f 

(bp)  (l+*  )  —  355 


z  U  3z2 


Combining  the  preceding  equations  with  equations  (55)  and  substituting 
the  resulting  equation  into  equation  (56)  give  the  following  principal 
equation  for  the  flow  on  a  surface  of  revolution: 


32\(f  (  wr2+wz2\  a2^  u  S\|r 

S^z  \1_  2  /  a  2+r  ^  M3z 


where 


V.lftf,  x  ,  Hz2,  ax 

az  az  a2  \  3z  r  dz 


N  =  -  (l+A  ) 


2s  fl  S  In  b  1  Ss*  .  1  1  Sll  a2-w2  f^X  1  / 1  Si  T  Ss 


cp  "rSp"  a2  r  5pf  a2  W'  w2\r^'r5^ 


For  this  equation,  the  characteristics'  are  real  or  imaginary  when  the 
resultant  relative  velocity  W  is  supersonic  or  subsonic,  respectively, 
For  the  subsonic  case,  it  is  again  better  to  use  the  following  form 
obtained  by  differentiating  equations  (55)  and  substituting  the  result¬ 
ing  equations  into  equation  (56) : 
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1  c)  In  bp  cty  9  In  bp 

2  <5cp  3q>+  9z 


27vnbp  - 


(bp)2  /  1  dl  T  dsN 
1  dt  V  "  r  ^ 
r^f> 


(59a) 


With  X  given,  and  h  determined  from  data  obtained  on  the  Sg 
surface,  there  are  now  the  six  independent  relations  equations  (21), 
(59)  or  (59a),  (55a),  (55b),  (53),  and  (41a)  for  the  determination  of 
the  six  main  variables  in  \[r,  Wu,  Wz,  Wr,  s,  and  I.  The  f-components 

are  not  involved  in  the  calculation.  If  the  flow  is  adiabatic  with 
uniform  I  and  s,  the  equations  are  further  simplified. 


* 


Flow  along  Cylindrical  Surface 

If  the  flow  near  the  walls  of  an  axial-flow  turbomachine  can  be 
considered  to  take  place  on  a  cylindrical  surface,  then 


nu  =  nz  =  fu  «  f z  =  Wr  0  (60) 

Equation  (35)  now  gives 


c)Wr 

c  a  -  -r: - 

dr 

'  (61) 

which  is  relatively  small,  (if  c  is  negligible,  b  can  be  taken  as 

1  everywhere.)  For  flow  without  change  in  radial  distance,  the  quantity 

I  can  be  replaced  by  Hw  ^  =  h  +  The  equations  governing  the 

cylindrical  flow  are  then  (compare  reference  29) 

%  =  rbpWz 

(62a) 

3?  “  “  bpWu 

(62b) 

(63) 

D 

Dt 


=  T 


Ds 

Dt 


(64) 
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where 


1  52\|r 


0wuwzl52*  /  ^VV+NM  +  M^.n 
-2  —  r®+\1-a2;Sz2  +  r^+M3i  ^ 


,,  5  In  h  ,  5s*  1 

M - sr-  +  -‘;t  sr 

a 


_  ,  1  5  In  b  1  5s*  1 

r  5^>  r  5cp  „2 


a2-Wu2  !  5^  a2-W2  T  5s 


52\|r  52i|r 

L-L--  1  -J-  rr“ " 


1  5  In  bp  5\|r  .  5  In  bp  5i|r 
“2  5cp  3cp +  5z  3z 


(M2/: 

l  3f  V 

r  5cp 


T  5s 
r  3® 


(65a) 


In  general,  the  circumferential  derivatives  of  and  s  are  to 

be  determined  by  the  inlet  flow  and  equations  (63)  and  (64)..  For  adia¬ 
batic  flow  with  uniform  Hw  and  s  upstream  of  the  blade  row,  these 
derivatives  are  equal 'to  zero  everywhere,  making  the  problem  much  sim¬ 
pler.  The  main  difference  between  this  simplified  case  and  the  ordinary 
two -dimensional  flow  on  a  cylindrical  surface  is  the  inclusion  of  the 
factor  b  in  equations  (62)  and  (65)  (in  general,  b  is  a  function  of 
cp  and  z) .  If  the  velocity  diagram  is  such  that  there  is  considerable 
radial  gradient  in  the  radial  velocity  or  considerable  variation  of  the 
distance  between  the  adjacent  streamlines,  the  factor  b  is  not 
negligible. 


EQUATIONS  GOVERNING  FLUID  FLOW  ON  S1  SURFACE  FOR  RADIAL-FLOW 
AND  RADIAL-DISCHARGE  MIXED-FLOW  TURBOMACHINES 


Flow  along  General  Sp  Surface 

For  turbomachines  with  radial  discharge,  r  and  cp  are  considered 
as  the  two  independent  variables;  that  is, 

q  =  q[]r,  cp,  z(r,cp[] 
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Then 


8q  _  dq  dq 

8r  ~  c)r  -  nz 

1  8q  _  1  dq  nu  dq 
r  8cp  r  3<p  nz  7z 


(66) 


2a  -  w  2a  4.  2a 

Dt  r  9r  r  8cp  J 

Equations  of  continuity  and  motion.  -  By  the  use  of  these  relations, 
the  equations  of  continuity  and  motion  become 


1  3(pWrr)  1  3(pWu) 


8r 


3cp 


«  pc* 


(67) 


where 


1  (  ~Mr  ~"u 

HTV 


Swr  Bwu  Swz\ 
pr  3F-  +  Hu  ST  +  nz3T/ 


(68) 


and 


W 


wrwu 


.  3WU 

Wu8Wr 

"  W: 

8WZ 

u  97" 

+  T 

’  8cp 

z  8r 

8WU 

awr 

WZ 

8WZ 

..  -  .  .  1 

8r 

r 

8cp~ 

r 

acp 

DWZ 

■  =  wr 

aw 

z+^u.; 

Dt 

3r 

r  < 

81 


-  a*»u  «  -  ^  +  T  U  +  fr  (69a) 


8r 


8r 


1  81  T  8s 
r  8cp  +  r  8cp 


(69b) 

(69c) 


with 


f 


1  (bh  ds\  1  1  5p 


(70) 


Principal  equation.  -  If  a  variable  b  *  is  introduced  such  that 

D  In  b1 


Dt 


(71) 
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C  *  &X  S* 


Tf  to- 


(7  la) 


in  which  the  integration  is  performed  along  a  streamline  on  the  surface, 
then  the  continuity  equation  (67)  can  he  written 

3(b'pW  r)  3(h'pW„) 


and  a  stream  function  \|r  can  he  defined  on  the  surface  with 


H=  -b'pw, 


(73a) 


Si  =  rbtpw 
3cp  p  r 


(73b) 


Here  b'  can  he  interpreted  as  the  thickness  (in  the  z-direction)  of 
the  stream  sheet  whose  mean  surface  is  the  Sq  surface  considered.  The 
continuity  equation  (72)  can  again  he  obtained  by  equating  the  mass  flow 
into  and  out  of  an  element  of  the  stream  sheet  as  defined  by  two  axial 
planes  dcp  angle  apart  and  two  cylindrical  surfaces  dr  distance  apart 
as  shown  in  figure  7(b).  As  before,  the  difference  of  \jr  at  any  two 
points  on  the  Sq  surface  is  equal  to  the  mass  flow  across  any  line 
connecting  these  two  points.  By  the  use  of  the  preceding  two  equations 
and  the  relation 


h  -  I  + 


T-  "  *  <Vp> 


+  (!&. 

\r  acp, 


the  principal  equation  for  the  flow  of  this  surface  is  obtained  from 
equation  (39b): 


n  Wr2\  92^  WrWu  q  a2t  / 

.  "  n2  /  a_2  "  '  2  r  3r3Cp  +  V1 


/  r 


L  +  M  +  n  at 

,2  8Cf£  M  ar  +  r  3$ 


where 


„  a  in  b'  ,  as*  .  i  ai  .  „ 

M  - - =— -  +  r—  +  -77  -  7TT  +  W, 


awz  a'^-(wr2+Wu2)  GD2r2+Wr2 


3r  3r” 
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1 

r 


9  In  b'  .  1  9s*  ,  1  l  31  .  „  8wz\ 

acp  +  r  3CP  a2r  acp+  wz  acp  J 

a2-(wr2^u2)  /  !  9i  T  as 

a\2  \  r  8CP  r  9CP  u  " 


+ 


2o)Wr 


Wz  ^z\ 

r  acp/ 


This  equation  is  seen  to  he  hyperbolic  or  elliptic  when  V^r^u  is 
greater  or  less  than  the  speed  of  sound,  respectively.  For  the  elliptic 
case,  it  is  preferable  to  use  the  following  form: 


a2i(f  _±_  agjr  /a  in  b Tp  a\|r  i  s  in  b!p  af\ 
dr2  +  r2  acp2  ~  \  3r  ar  +  r2  acP  3<y 


(*>'p)2  r 

i21  L“ 

r  acp 


+  +  Wz  aWz 

r  3cp  +  r  3cp  +  u  +  r  acp 


=  0 


(75a) 


The  integrability  condition  (51)  is  now  written 


hence 


(76) 


(76a) 


The  procedure  of  solving  the  principal  equation  with  the  various  terms 
in  it  determined  by  other  flow  equations  is  the  same  as  that  in  the 
previous  system. 


Flow  along  a  Surface  of  Revolution 


For  the  special  case  of  flow  on  a  surface  of  revolution,  equa¬ 
tions  (52)  and  (53)  hold  (with  X  considered  as  a  function  of  r)  and 
the  expression  of  c!  reduces  to 


c  1 


1  dwr  _  dwz 

X  <3z  ”  dz 


(77) 


POAn 
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Flow  on  Radial  Plane 

For  the  special  case  of  flow  on  a  radial  plane, 

Hr  “  nu  =  fr  =  fu  “  wz  “  0 
and  equations  (77)  to  (81)  reduce  to 

, 

c  “  ■  sr 


u 


h  -  I  + 


l^.rb'pWr 


1  2  2 

2“r 


(82) 

(83) 

(84a) 

(84b) 

(85) 


=  0 


(86) 


where 


N  =  -  - 


M  = 

S  In  b' 

L  Ss*  1 

dr 

'  +  ^  +  ?" 

1  s_ 

In  b' 

1  ds* 

+  1  /a2-"u 

r 

"Sep 

r  W 

a2r  Wr2 

.2-W2-Wr2-cu2r2 


Wi 


a2  Wr 


and 


d^|r  1  1  /S  In  b'p  1  d  In  b 1  p  cty  \  ,  ; 

5?  +  ?  *  +  ■?  ^  ®  +  7  “  p 


(b'p) 
l  5\j? 


(i  Si  t 
\r  3<p"r 


Ss 

^}m  ° 


(86a) 
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Alternative  Form  of  Equations  for  Flow  along 
Surface  of  Revolution 

i 

The  equations  given  in  the  preceding  sections  are  obtained  for 
turbomachines  to  avoid  an  infinite  value  of  the  partial  derivative  with 
respect  to  z.  Difficulty  still  exists  in  using  either  of  the  systems 
in  the  case  of  a  mixed-flow  type  machine  with  an  axial  inflow  and  a 
radial  outflow.  For  solutions  of  general  Sq  surfaces,  this  difficulty 
can  be  avoided  by  dividing  the  machine  at  the  middle  of  the  flow  path 
and  using  the  first  system  at  the  inlet  portion  and  the  second  system 
at  the  exit  portion.  If  the  Sq  surface  can  be  approximated  by  a  sur¬ 
face  of  revolution,  it  is  convenient  to  use  a  set  of  orthogonal  coordi¬ 
nates  l  and  cp,  where  l  is  the  arc  length  of  the  generating  line  of 
the  surface  of  revolution  in  the  meridional  plane  and  Cp  is  the  usual 
cylindrical  angle  (fig.  8) .  Because 

^r  dr  , n  N 

—  s  —  s  sm  or  (87} 


WZ  ~  d l 


-  cos  a 


then,  for  use  with  the  first  system, 


8q  5q 

Esseca5? 


and,  for  use  with  the  second  system. 


3q  3q 

—  »  CSC  O 

3r  dZ 


By  use  of  the  preceding  relations,  the  equation  of  motion  in  the  circum¬ 
ferential  direction  as  given  by  either  equation  (56)  or  (79)  for  the 
two  systems,  respectively,  becomes  in  both  cases 


1  owZ 
r  3V 


-l  -+  2m 


sm  a 


1  hi 
r  <5cp 


T  Ss 
r  dCD 


(79a) 


which  agrees  with  the  results  obtained  in  references  29  and  30  in  a 
different  manner.  The  subsequent  equations  given  in  these  two  refer¬ 
ences  can  be  modified  and  used  for  such  surfaces.  (The  last  term  on  the 
left  side  of  equation  (79a)  represents  the  rotationality  of  the  absolute 
flow  and  is  not  included  in  reference  29,  which  is  derived  for  irrota- 
tional  absolute  flow.) 
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By  comparing  the  integrating  factor  h  used  herein  and  the  thick¬ 
ness  of  the  stream  filament  of  revolution  t  used  in  these  two  refer¬ 
ences,  it  is  seen  that  the  two  play  exactly  the  same  role  in  the  con¬ 
tinuity  relation.  Although  b  is  obtained  mathematically  as  an  inte¬ 
grating  factor,  physically  it  may  then  be  visualized  as  the  thickness  - 
of  the  stream  filament  in  the  r  or  z-direction  for  the  two  systems , 
respectively.  The  use  of  b  herein  is,  of  course,  more  general  in  that 
it  varies  two -dimensionally  over  the  surface  in  the  general  case,  where, 
as  in  references  29  and  30,  t  is  considered  a  function  of  l  only. 


EQUATIONS  GOVERNING  FLUID  FLOW  ON  S2  SURFACES 


In  the  preceding  section,  it  was  shown  that  the  determination  of 
the  flow  on  Sq  surfaces  requires  a  knowledge  of  the  radial  variation 
of  the  velocity  components.  This  knowledge  can  be  obtained  by  following 
the  fluid  motion  along  relative  stream  surfaces  of  the  second  kind,  S2. 
On  S2,  the  relations  (28)  to  (31)  also  hold.  These  relations,  how¬ 
ever,  will  now  be  used  to  eliminate  the  independent  variable  cpj  that 
is,  any  quantity  q  on  S2  is  now  considered  as 

q  =  q[^r,z,cp(r,z)] 

Accordingly,  on  S2 


3q  _  Sq  \  1  5q 
3r  “Sr  nu  r  Sep 

3q  _  Sq  nz  1  Sq 
3z  “  Sz  "  nu  r  Sep 


(91) 


Equations  of  Continuity  and  Motion 


Equations  (30)  and  (91)  are  used  to  change  the  continuity  equa¬ 
tion  (35)  to 


1  a(pwrr)  3(pWz) 
r  3r  +  8z 


P  C(r,z) 


(93) 


* 
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where 


3W, 


aw, 


aw. 


/  \  *1  (  WH17 

°(r,z)  ■"  "  —  i  nr  3^-  +  nvi  5^-  +  nz  ^cp- 


nur 


(94) 


For  general  rotational  motion,  the  equations  of  motion  (14a)  in  the 
three  perpendicular  directions  are 


W  ( 

_ 
r  L" 


d(Vur) 

^r 

/awr 

_  dr 

■  WJ  z 

\3r  ■ 

W 


Pb(Vur)  awrl  / 

L~ ■  wj  " Wz  v 

aw, 


-  w„  i- 


/aw. 


r  \3z  dr  y 


9- 

i)..i 


a?  +  T  3F 


3s  ~\ 


aw 

>Cf 

aw 


1  owz  awu' 

r  acp  az 


ai  ■  t  as 

<5cp  +  r  3cp 

as 


>  (95) 


33  +  T  33  J 


In  following  the  motion  on  Sgj  equations  (95)  are  reduced  to  the  fol¬ 
lowing  form  by  using  equations  (9),  (16),  (31),  and  (91): 


^  a(  V) 

r 


/  awr  aw_\ 

3r  +Vz^--—)  =  - 


|£  +  T  +  Fr 
3r  ar  r 


wr  9(vur)  wz  a(vur) 

r  9r  +  r  az  “ 


u 


or  Fur  ss 


D(Vur)l 
Dt  J 


-  W,. 


/  3Wr  9WZ\  Wu  9(Vur) 


*  \az 


ar 


az 


+  T  —  +  F 
az  az  +  z 


(96a) 

(96b) 

(96c) 


where  F  is  a  vector  having  the  unit  of  force  per  unit  mass  of  gas 
defined  by: 


b  — —  (d  -  t  dV  *  — —  - 

nur  \dcp  acp/  nur  p 


aP 

35  n 


(97) 


A  similar  result  is  obtained  for  the  equation  of  motion  in  the  form  of 

(2): 

2 


TT  .  tt  3Wr  VU  _  1  3p 

r  3r  2  3z  "  r  “  '  p  ar  Fr 


a 


»r  a(vur)  wz  a(yuri 


r  ar 

aw, 


+  — 
r 


cjw7  aW7  1  am 

wz  — £  +  wr  —5.  -  -  i  §E  +  Fz 
z  3z  r  ar  p  3z  z 


az 


u 


j 


(98) 
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Because  the  vector  F  is  normal  to  the  S2  surface, 

FrWr  +  FUWU  +  FZWZ  =  0  (99) 

By  the  use  of  equations  (99)  and  (96),  it  can  be  shown  for  steady  flow 
on  an  S2  surface  that 


DI  _  T  Ds 
Dt  Dt 


(99a) 


This  result  is  the  same  as  that  obtained  for  the  Sq  surface.  Again, 
for  the  present  problem  of  steady  relative  flow  on  the  S2  surface,  the 
relation  (99a)  can  be  taken  either  as  one  of  the  equations  of  motion  or 
to  represent  the  relation  (99).  In  other  words,  there  are  only  four 
independent  relations  among  equations  (96a),  (96b),  (96c),  (99),  and 
(99a) .  In  the  following  development,  it  is  found  convenient  to  use 
equations  of  motion  in  the  form  of  equation  (96),  not  only  because 
bl/br  is  zero  in  many  design  problems  (whereas  dp/dr  ^  0) ,  but  also 
because  equation  (96)  leads  to  a  form  capable  of  a  rigorous  solution 
for  both  subsonic  and  supersonic  flow  and  shows  clearly  how  the  various 
design  factors  affect  the  three-dimensional  motion  in  general.  (See 
equations  (106)  to  (114)  that  follow.) 


In  a  manner  analogous  to  the  Sq  surface,  the  continuity  equa¬ 
tion  (93)  is  put  into  the  form 


3(rBpWr) 

Sir 


8(xBpWj 


(100) 


by  the  use  of  an  integrating  factor  B,  which  is  related  to  C  by  the 
following  equation; 


or 


D  In  B  TT  3  In  B  ,  TT  3  In  B 
“Dt  -  =  Wr  “ S»r  +  -dT~ 


(101) 


(101a) 


Equation  (100)  is  the  necessary  and  sufficient  condition  that  a 
stream  function  \|r  exist  and 


U  -  rBpWz 


(102a) 
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If  -  -  rBpW r  (102b) 


The  difference  in  ¥  at  two  points  j  and  k  on  the  Sg  surface  is 

fk  Pk 

¥  -  ¥J  » J  &f  =  I  rBp(Wz  dr  -  Wr  dz) 


Similar  to”  the  flow  on  the  Sq  surface,  the  preceding  equation 
indicates  that  B  is  proportional  to  the  angular  thickness  of  a  thin 
stream  sheet  whose  mean  surface  is  the  stream  surface  S2  considered 
herein  and  whose  variable  circumferential  thickness  is  equal  to  rB. 
Indeed,  if  the  mass  flow  into  and  out  of  the  element  of  such  a  stream 
sheet  (cut  between  two  planes  normal  to  the  z-axis,  and  a  distance  dz 
apart  and  between  two  cylindrical  surfaces  dr  apart  (fig.  7(c)))  is 
equated  to  zero  and  the  distances  dr  and  dz  approach  zero  as  a  limit, 
the  following  equation  is  obtained: 


a(TpWr)  3(tP¥z) 
8r  +  dz 


(100a) 


Comparing  this  equation  with  equation  (100)  and  considering  the  mass 
flow  relations  show  T  to  be  proportional  to  rB.  This  proportionality- 
means  that  B  can  be  physically  intepreted  as  a  quantity  which  is  pro- 
portional  to  the  angular  thickness  of  a  stream  sheet  whose  mean  surface 
is  the  Sg  surface  considered  herein.  With  this  interpretation,  B  is 
immediately  seen  to  be  closely  related  to  the  angular  distance  between 
two  neighboring  blades.  In  actual  calculation,  only  the  ratio  rB  to 
(rB)q  or  T  to  T-[  is  important,  and  it  is  also  easier  to  obtain  the 
variation  in  rB  'from  the  distance  between  adjacent  streamlines 
obtained  on  Sq  surfaces  than  to  evaluate  B/Bq  by  equations  (101a) 
and  (94)  using  data  obtained  on  Sq  surfaces. 


Principal  Equation  for  Case  with  Vur  Given 

In  the  solution  of  flow  on  an  S2  surface,  the  continuity  equations 
and  the  equation  of  motion  in  the  radial  direction  are  combined  to  form 
the  principal  equation.  The  principal  equation  will  now  be  obtained  for 
two  main  groups  of  present  designs  ’in  which  a  certain  desirable  varia¬ 
tion  of  the  angular  momentum  of  the  fluid  Vur  and  of  the  ratio  of 
relative  tangential  and  axial  velocity  are  prescribed  on  the  S2  m  sur¬ 
face,  respectively.  These  equations  can  also  be  used  for  the  solution 
of  a  direct  problem,  in  which  the  same  information  obtained  on  Sq 
solutions  of  a  previous  cycle  is  used  as  known  values  in  the  S2 
solution. 
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For  the  first  group,  the  following  equation  is  considered  known: 

Vur  =  G(r,z)  (103) 

Among  this  group  of  designs  are  the  free -vortex  design  (in  which  G  is 
simply  a  function  of  z),  the  more  general  "solid-body  rotation"  design, 
the  "symmetrical  velocity  diagram  at  all  radii"  design,  and  others  (foi 
example,  see  references  17  and  18) . 


From  equations  (102)  and  (45), 

3W.  a2i 


rBp 


f_!  _  £_I  /  1  _1_ 

ar  ~  2  +  l "  r  ~  2 

8r  \  a 


ah  as^  a  in  b\  aw 
3r  +  3r  ”  3r 


i\  aw 

/  ar 


rBp 


aw. 

az" 


-  ii  /  _i_  3h  as^_  a  In  b\  aw 

”  ^_2  +  V  .2  3z  +  3z  3z  )  3z 


3z  \  a 
But  from  equations  (9)  and  (101), 

T  ,  o>2r2  Wu2  1  ,  ^  _2 

h  =  I  +  — ^ - o - o’  (rBp) 


(104a) 

(l04h) 


2  2  2 
Differentiating  with  respect  to  r  and  z  gives 


(105) 


a2-(W  2+W  2) 
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ar  ar 
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(rBp)-l^Wz  -  Wr  |^) 


Substituting  the  preceding  equations  into  equations  (104)  and  adding 
give 


[a2-  (wrW)]  -BP  (£  -  £)  =  (aB-WrB) 
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Substituting  equation  (106)  into  equation  (96a)  and  dividing  by 
a2  yield  the  following  principal  equation  for  the  fluid  flow  on  sur¬ 
face  Sg: 


a^  /  8r^  a^  3r3z  \ 


a2w 

3z2 


3¥ 

+  'N  —  + 

3r 


3¥ 

M  =  0 

3z 


(107) 


where 


3  In  B  3s^_  _  _1_  ( 31  _  w  3Wu\ 

3z  3z  “  a2  \3z  u  3z  / 
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1 
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3  In  B  3s*  1  /  31 

3r  +  ar  “  a  2  \3r 


"  wu  +  “2r 


+ 


a2-(wr2+Wz2) 

•  aV. 


m  3S  _ 
1  3r  + 


Wu  3(Vur)" 
+  r  3r  _ 


From  the  coefficients  of  the  second  derivatives,  the  principal  equation 
is  seen  to  be  hyperbolic  or  elliptic  when'  the  meridional  velocity 

W j  «  /\Jwr2  +  Wz2  is  greater  or  less  than  the  speed  of  sound,  respec¬ 
tively,  For  the  elliptic  case,  it  is  again  convenient  to  write  the 
principal  equation  in  a  slightly  different  form.  From  equation  (101), 

rCp  3Wz  a2y  1  3W  3  In  Bp  g¥ 

H  3r  ar2  r  3r  3r  3r 

8Wr  91n  Bp 

-  rBP  =  — 2  -  "  a  z  ' 

3z 

Substituting  into  equation  (96a)  results  in 
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(107a) 
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With  the  variation  of  Vu  or  Wu  prescribed  by  the  designer  in  an 
inverse  problem  or  taken  from  the  previous  calculation  in  a  direct 

problem,  the  meridional  velocity  components  are  determined  by  equa¬ 
tions  (107)  and  (l07a) .  (Other  equations  are  used  to  determine  various 
terms  involved  in  the  coefficients  M  and  N.) 


Principal  Equation  for  Case  with  W^W,^  Given 

In  the  second  group  of  designs,  the  following  relation  is  pre¬ 
scribed  on  an  surface  (for  example,  see  references  17  and  18): 


(109) 

In  order  to  result  in  blades  with  the  mean  blade  surface  composed  of  all 
radial  elements  (for  high-speed  rotation),  it  maybe  desirable  to  spec¬ 
ify  a  mean  S£  surface  consisting  of  all  radial  elements.  Then 

ir=r6l(z)  (110) 

Similarly,  in  order  to  obtain  a  cooled  turbine  rotor  blade  with  minimum 
twist,  the  following  function  may  be  specified  on  S2^m: 

Wu 

%  -  e2(z>  <m> 


In  application  to  direct  problems,  one  of  the  preceding  relations  is 
obtained  from  the  solution  in  the  previous  cycle  and  is  considered 

as  given  in  the  S2  solution.  In  both  inverse  and  direct  problems, 
with  the  relation  between  Wu  and  Wz  given  by  these  equations,  all 

three  velocity  components  are  to  be  combined  into  the  main  terms  of  the 
principal  equation  as  follows:  Substituting  relation  (109)  into  equa¬ 
tion  (96a)  gives 


<1+s2)  ^  +  g  (f  +  g )  W,  +  a»s  +  i  {-  g  +  T  fi  +  Fr)  -  0 


(112) 


Instead  of  equation  (105) , 
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should  now  he  written.  Differentiating  with  respect  to  r  and  z, 
combining  with  equation  (104),  and  substituting  into  equation  (112) 
give  the  following  form  of  the  principal  equation: 


<^>  VS)  &  -2<i+b2) 


WrWz  d2y  (  Wu2+Wz2\  a2^r  '  8f 


2  3r3z 


where 


+  VI- 
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3r  3z 
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-+i( 
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3r  z  &  Sr 


si"!  g  /s  +  |s\  +  ^  (.  |i  +  T  |i  +  pr  - 
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This  equation  is  hyperbolic  when  the  relative  velocity  is  supersonic, 
elliptic  when  the  relative  velocity  is  subsonic.  For  the  subsonic  case, 
a  form  of  this  equation  more  convenient ' for  computation  is  obtained  by 
substituting  equation  (108)  into  (112): 


3g\  3¥  aff 

ai7  8r  8z2  ’ 


[(1+S3)1^|I+^|I] 


2@corBp  +  hgl!  (-  ||  +  T  ff  +  Fr^=  0  (114a) 

3r 

It  may  be  noted  that  for  both  groups,  equation  (96a)  rather 
than  (96c)  is  chosen  to  obtain  the  principal  equation  of  the  present 
problem,  because  Fr  is  always  much  smaller  than  Fz  in  axial 
machines  and  Fr  is  zero  or  nearly  zero  on  Sg  surfaces  for  high-speed 
centrifugal  and  mixed-flow  impellers  whose  mean  blade  surfaces  are 
usually  composed  of  all  radial  elements.  (For  low-speed  centrifugal 
impellers,  equation  (96c)  can  be  used  to  form  the  principal  equation 
in  a  similar  manner.) 
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Procedure  of  Solution 

Although  the  equation  of  motion  (96a)  is  chosen  to  form  the  prin¬ 
cipal  equation,  other  equations  are  to  he  used  to  obtain  the  various 
terms  involved  in  the  principal  equation.  As  in  the  case  of  general 
Sq  surfaces,  there  are  ten  basic  variables  to  define  the  flow  and  the 
shape  of  the  Sg  surface.  They  are:  y,  B,  Wr,  Wu,  Wz,  Fr,  Fu,  Fz,  s, 
and  I  (or  p).  B  is  considered  given,  (in  the  direct  problem,  B 
is  evaluated  directly  from  the  distances  between  adjacent  streamlines 
or  according  to  equation  (100a)  using  the  value  of  C  obtained  on  Sq 
surfaces;  in  the  inverse  or  design  problem,  B  is  estimated  (refer¬ 
ences  29  and  35)  from  the  blade  thickness  as  desired  from  blade  stress 
and  other  considerations.)  On  the  other  hand,  there  are  seven  inde¬ 
pendent  relations  in  one  energy  equation  (21);  three  equations  of 
motion,  one  of  equations  (107),  (107a),  (114),  or  (114a),  and  equations 
(96b)  and  (96c);  the  orthogonality  relation  between  W  and  F,  equa¬ 
tion  (?9a);  and  the  two  equations  relating  V  and  velocity,  (102a) 
and  (102b). 

Direct  problem.  -  In  the  direct  problem,  two  alternative  procedures 
may  be  used.  If  the  shape  of  the  S2  surface  (determined  from  the 
data  obtained  on  Sq  surfaces)  is  considered  as  given  in  the  present 
S2  solution,  two  additional  relations  between  the  n-  or  F-components 
completely  define  the  problem.  The  procedure  of  calculation  is  as 
follows: 

(1)  Use  equations  (20)  and  (21)  to  determine  the  variation  of 
s  and  I . 

(2)  Compute  Wu  from  the  orthogonality  relation  as  follows: 

wu  =  -  (  —  Wr  +  —  w  \ 

u  \nu  r  %  7 

(3)  Compute  Fu  from  equation  (96b). 

(4)  Solve  the  principal  equation. 

(5)  Obtain  Wr  and  Wz  from  equations  (102). 

If  only  the  tangential  velocity  or  the  relation  (109)  is  taken 
from  the  Sq  solutions  of  the  previous  cycle  and  is  considered  as 
given  in  the  present  Sg  solution,  one  more  relation  is  available  * 

between  the  F-components  such  as  that  which  exists  between  the 
f-components  on  the  Sq  surface: 


F  •  7  X  F  =  0 


(115) 
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Writing  equation  (115)  in  scalar  form  and  using  the  relations  (3l) 
and  (91)  give 


By  integrating  along  a  constant  z-line,  equation  (115a)  provides  the 
following  relation  to  determine  the  value  of  Fr  to  he  used  in  the 
principal  equation  from  the  values  of  Fu  and  Fz: 


■^uVo 


8.  /Fz  \  . 

8r  \FurJ  dz 


(115b) 


If  Fr  =  0,  at  zQ 


Fr  =  Fur 


(?z_\ 

WV 


(115c) 


The  procedure  of  calculation  is  as  follows: 


(1)  Use  equations  (21)  and  (99a)  to  determine  the  variation  of  s 
and  I. 

(2)  Compute  Fu  and  Fz  from  equations  (96b)  and  (96c). 

(3)  Compute  Fr  from  equation  (ll5b)  or  (115c). 

(4)  Solve  ¥  from  the  principal  equation. 

(5)  Compute  Wr  and  Wz  from  equations  (l02a)  and  (l02b). 

Inverse  problem.  -  In  the  inverse  or  design  problem  of  a  finite 
number  of  thick  blades,  in  addition  to  the  blade-thickness  distribution 
or  its  equivalent  B,  either  equation  ( 103)  or  (109)  is  prescribed  on 
a  mean  stream  surface  It  may  appear  that  still  another  rela¬ 

tion  can  be  prescribed  on  the  mean  surface.  The  differentials  of  the 
coordinates  of  the  surface  are  now  governed  by 


Fr  dr  +  Fur  dcp  +  Fz  dz  =  0 


(116) 


48 


NACA  TS  2604 


however,  and,  in  order  that  this  differential  equation  will  lead  to  an 
integral  surface  of  the  form  represented  by  equation  (28),  F  must 
satisfy  the  condition  of  integrability  as  given  by  equation  (115) 
(reference  34).  An  expression  similar  to  equation  (ll5a)  for  the  case 
of  an  infinite  number  of  blades  was  first  pointed  out  by  Bauersfeld 
(reference  2)  in  a  discussion  of  the  Lorenz  paper  (reference  l) .  In 
effect,  it  restricts  the  freedom  that  the  designer  has  in  prescribing 
the  velocity  components  of  the  fluid  on  the  surface.  Hence,  in  the 
inverse  problem  of  a  finite  number  of  thick  blades,  in  addition  to  the 
blade  thickness  distribution  or  its  equivalent  B,  the  designer  can 
specify  only  one  relation  on  the  mean  stream  surface,  which  relation  may 
be  either  the  tangential  velocity  as  given  by  equation  (103),  the  flow 
angle  between  the  tangential  and  axial  velocity  as  given  by  equa¬ 
tion  (l09 ) ,  the  axial  velocity,  or  any  other  reasonable  relation  that 
will  lead  to  a  solution  of  the  set  of  equations. 

In  the  preceding  consideration,  the  hub  and  casing  shapes  are  also 
prescribed  by  the  designer  in  the  inverse  problem.  Alternatively,  the 
prescription  of  the  hub  shape  can  be  replaced  by  a  prescription  of 
another  relation  at  the  casing,  thereby  fixing  the  shape  of  and  the  flow 
along  S2,m  at  the  casing  entirely.  The  flow  is  then  extended  to  the 
hub  and  the  last  streamline  gives  the  hub  contour  (reference  19). 


Approximations  Involved  in  Through-Flow  Theory 

When  the  equations  previously  derived  in  reference  18  for  a  large 
number  of  thin  blades  are  compared  with  the  corresponding  equations 
derived  herein  along  a  stream  surface,  the  two  are  obviously  exactly 
the  same  if  the  ordinary  derivatives  used  in  reference  18  are  replaced 
by  the  present  partial  derivatives  following  the  stream  surface,  and 
if  B  is  equal  to  1  or  if  the  variation  of  B  along  the  flow  path 
is  zero.  In  the  interpretation  of  the  through-flow  solutions  as  the 
flow  along  a  mean  stream  surface  (which  divides  mass  flow  into  two  equal 
parts  circumferentially)  or  as  the  flow  along  the  mean  channel  surface 
(geometrical  mean),  the  first  difference  can  easily  be  removed  by  simply 
interpreting  the  values  obtained  in  the  solution  as  those  along  the  sur¬ 
face  rather  than  in  the  meridional  plane.  The  second  condition,  however, 
is  satisfied  only  when  the  circumferential  variation  of  all  the  velocity 
components  approaches  zero,  or  when  the  circumferential  derivative  of 
the  tangential  velocity  and  the  ratios  of  nr  and  nz  to  nu  approach 
zero  (see  equation  (94)). 

Besides  the  use  as  a  limiting  solution  in  general  and  to  give  cer¬ 
tain  trends  where  the  contribution  due  to  the  finite  number  of  blades 
is  small  or  constant,  the  through-flow  calculation  should  be  properly 
modified  by  the  factor  B  in  its  application  to  actual  turbomachines 


MCA  IN  2604 


49 


of  a  finite  number  of  thick  blades.  As  B  can  he  physically  inter¬ 
preted  as  the  ratio  of  the  local  angular  thickness  of  the  stream  sheet 
to  its  inlet  value,  a  good  approximate  value  can  he  obtained  by  solving 
the  two-dimensional  flows  on  a  number  of  stream  surfaces  of  revolution 
starting  at  different  radii  at  the  inlet.  For. the  subsonic  flow  in  the 
turbine  cascade  reported  in  reference  29  and  for  the  supersonic  flow  in 
two  impulse  bladings  investigated  in  reference  30,  the  reductions  in 
angular  thickness  from  the  inlet  value  along  the  mean  streamline  are 
seen  to  be  a  chordwise  average  of  4  and  9  percent  more  than  the  reduc¬ 
tion  in  the  channel  width,  respectively.  Also,  in  the  subsonic  case, 
the  influence  is  extended  a  certain  distance  outside  the  blade  row.  The 
inclusion  of  this  factor  B,  even  if  it  is  approximate,  should  give  a 
much  closer  answer  than  that  obtained  with  B  taken  as  1. 

In  this  interpretation  of  the  infinite  number  of  blades  solution 
as  the  solution  of  through  flow  along  a  particular  stream  surface 
between  two  adjacent  blades,  the  distributed  "body  force"  F  has  a 
definite  meaning,  as  given  by  equation  (97).  (For  an  infinite  number 
of  blades,  F  becomes  the  blade  force.)  For  blades  with  large  turning 
and  large  radial  twist,  as  in  a  free-vortex  turbine,  the  influence  of 
the  radial  component  of  F  on  the  flow  is  not  negligible. 

CIRCXJMFEEEOTIAL  VARIATION  OF  FLUID  PROPERTIES  BY  USE  OF  POWER  SERIES 

In  general,  the  blade -to -blade  variations  of  fluid  properties  are 
to  be  obtained  from  calculations  on  Sq  surfaces.  When  the  twist  of 
the  Sq  surface  is  large,  some  other  method  of  obtaining  the  blade-to- 
blade  information  is  desirable.  For  subsonic  irrotational  absolute 
flow,  this  information  can  be  obtained  by  extending  the  solution 
obtained  on  the  mean  stream  surface  in  the  circumferential  direction  by 
the  use  of  power  series  (without  the  consideration  of  the  shape  of  the 
Sq  flow  surfaces).  The  various  derivatives  involved  in  the  series  are 
bbtained  from  the  flow  condition  on  the  mean  stream  surface.  The  higher 
the  solidity  and  the  thinner  the  blade  sections,  the  fewer  are  the  terms 
required  for  a  given  accuracy.  Results  obtained  in  references  29  and  36 
indicate  that  only  three  terms  in  the  series  will  be  required  to  give 
sufficient  accuracy  for  high-solidity  turbines  and  centrifugal  com¬ 
pressors. 

The  series  method  will  also  be  used  in  one  of  the  two  methods  of 
the  inverse  solution  in  which  the  flow  obtained  on  the  mean  stream  sur¬ 
face  is  extended  out  circumferentially. 

Denoting  the  absolute  vorticity  VXV  by  |  and  using  the  rela¬ 
tions  (16),  (9l),  and  (97)  in  equations  (l9)  give 
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1  toz  S(Vur) 
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From  the  preceding  equations, 
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aw 
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Fz  3(Vur)  Fr  3(Vur) 
Fur  3r  +  Fur  8z 


(118) 


This  equation  means  that  the  apparent  vorticity,  which  is  obtained  by 
differentiating  the  velocity  on  the  mean  stream  surface  with  respect  to 
the  coordinates,  is  not  zero  even  if  the  absolute  vorticity  is  zero  or 
tangent  to  the  mean  stream  surface.  Substituting  equations  (U7a) 
and  (117c)  into  equation  (94)  results  in 
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3(Vur)  Fz  3(Vur) 
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ar  Fur  az 

Substituting  equation  (119)  into  equations  (ll7a)  and  (ll7c)  gives 


1  Swz  =  1 

3(Vur) 
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Fur  3z 
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F_  Fr 


(121) 


The  first  derivatives  of  h  or  p,  and  p  can  be  obtained  as  follows 
From  equation  (97), 


1  dh  _  T  as  _  J-__  OT  _  „  p 
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1  d] 


(122) 
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Or  from  equation  (9), 

1  Sh 

r  Sep 


isi  /ira»  wu  aw  w2swz\  . 

=  r59"\T^_  +  T^'  +  7^7  (122a) 


With  ^  known,  can  he  obtained  by  using  equation  (12b): 


1  5  In  p  =  J_  Sh  _  1  Ss* 
r  a2r  5q5  r  Sep 


(123) 


The  second  derivatives  of  the  fluid  properties  with  respect  to  cp 
can  be  obtained  in  a  similar  manner.  Differentiating  the  continuity 
equation  (l)  with  respect  to  cp  and  dividing  by  r  give 


1  S2(pWrr)  ±  S2(pWu)  ]_  S2(pWz) 
~^2,  Sr  Sep  +  r2  +  r  Sz  Sep 

Equations  (9l)  are  used  to  change  equation  (124)  to 


(124) 
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1 

a  S(pWrr)  ^  Fr  S2(pWr) 

_  1 

a  S(Pwz)  (  fz  S2(Pwz) 

r2  Sep2 

"  r2 

_ar  Sep  ‘  Fu  ^cp2  _ 

r 

_8z  Sep  '  Fur  ^  _ 

(124a) 

Differentiating  equations  (ll7a)  and  (117c)  with  respect  to  cp  and 
dividing  by  r  result  in 

_1_  =  1  Hr  1_  5^Vur)  =  1  S^r  Q_  d(Vur)  ,1_  £z 

r2  ^  r  cSep"  +  r2  Sz  Sep  r  ScfT  r2  dz  Sep  r2  Fu 

(125) 

1  _  lHz  1  ^(Vur>  1  ^z  1  3  My)  l  Fr 

r2  ^2  "  r  Sep”  +  r2  Sr  Sep  ”  r  Sep  +  r2  3r  Sep  +  r2  Fu  ^,2 


(126) 
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Substituting  equations  (125)  and  (126)  into  (124a)  and  noting  that  F 
is  perpendicular  to  W  give 


1  d2wu 
r2  Sep2 


Fu2  f  Fr  _3  5(Vur)  Fz  5(Vur)  1  Q_  5(PWrr) 
F2  Fur2  3r  ^  +  V  3z  +  pr2  3r  ' 


1  a  d(pWz)  2C  S  In  p  Fz  d£r  Fr  d|2 
pr  az  cKp  r  5q>  +  Fur  <3cp  Fur  Sep  _ 


(127) 


Equation  (103)  is  to  "be  used  in  equations  (lOl)  and  (102)  to  obtain  the 
second  derivatives  of  Wr  and  Wz.  The  second  derivatives  of  h  and 

p  are  obtained  from  equations  (9)  and  (123)  as 


c^h  _  d£l  __  _1_  S2W2 
r2  Scp^  r2  Scp^  2r2  dcp£ 


(128) 


1  S2  In  p  1  a2  h  1  S2s  . 

r2  "a2r2Scp2  “  r2  Sq>2  U 

Similar  formulas  can  be  obtained  for  higher-order  derivatives.  At 
a  fixed  value  of  r  and  z,  the  velocity  components,  h,  and  p  at  a 
short  angular  distance  away  from  the  mean  stream  surface  can  then 

be  obtained  by  a  Taylor  series: 


(cp-qp  )2 

q(«P)  =  qOPm)  +  (<?>-%)  q'(<%)  +  — g2-  q'^cpm) 


+ 


*3.'  (^m)  +  •  •  • 


(130) 

An  alternative  way  to  obtain  density  is  to  use  equation  (145)  (to  be 
given  subsequently)  after  the  other  fluid  properties  are  determined. 
Obviously,  the  preceding  equations  are  most  useful  when  the  flow  is 
isentropic  with  vorticity  equal  to  zero.  Otherwise,  the  variation  of 
vorticity  along  the  mean  stream  surface  has  to  be  determined  first. 

At  present  no  such  method  is  available.  It  appears,  however,  that  the 
method  of  Squire  and  Winter  (reference  37)  and  Hawthorne  (reference  38) 
may  be  generalized  to  compressible  flow  for  the  variation  of  vorticity 
along  a  mean  stream  surface  in  turbomachines. 
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STEPS  FOR  COMPLETE  SOLUTIONS  OF  THREE-DIMENSIONAL 
DIRECT  AND  INVERSE  PROBLEMS 

In  general,  the  solution  of  the  three-dimensional  direct  and  inverse 
problem  involves  the  use  of  both  Q±  and  Sg  surfaces.  In  the  direct 
problem,  starting  with  assumed  flow  surface,  the  solution  is  obtained 
through  the  successive  (alternate)  use  of  the  two  kinds  of  flow  sur¬ 
face,  although  a  satisfactory  approximate  solution  may  be  obtained  in 
one  or  two  complete  cycles.  The  use  of  an  approximate  method  of  solu¬ 
tion  to  get  a  good  starting  value  on  each  surface  will  shorten  the  length 
of  computation.  For  inverse  problems,  the  process  is  usually  shorter. 

The  calculation  will  start  on  the  S2^m  surface  on  which  either  a  con¬ 
dition  on  the  fluid  velocity  or  the  shape  itself  is  prescribed  and  an 
estimated  value  of  B  for  a  desirable  blade  thickness  distribution  is 
used.  After  the  solution  on  the  S2,m  surface  and  its  shape  are  obtained, 
the  blade  coordinates  are  obtained  by  extending  the  solution  circumferen¬ 
tially  either  by  the  series  method  or  by  the  method  given  in  reference  35 
using  the  variation  of  the  distances  between  the  streamline  obtained  in  the 
£>2,m  surface.  Because  it  is  important  only  to  obtain  the  right  order  of 
magnitude  and  the  right  kind  of  variation  (three -dimensionally)  of  the 
blade  thickness,  the  first  solution  may  give  satisfactory  results.  The 
velocity  distribution  on  the  blade  surface  is  controlled  directly  by  the 

one  relation  specified  on  the  Sg  m  surface  and  the  variation  of  B. 

« 

Suitable  procedure  is  subsequently  suggested  for  the  solutions  of 
direct  and  inverse  problems  with  either  irrotational  or  rotational  inlet 
absolute  motion,  at  design  or  off -design  flow  conditions,  for  turbo- 
machines  having  various  wall  configurations  (fig.  6). 


Direct  Problem 

Axial  turbomachines  with  nontapered  straight  walls.  -  In  this  type 
of  machine,  it  is  desirable  to  start  the  computation  on  Sp  surfaces, 
because  with  short  axial  blade  length,  the  total  deviation  of  the  Sq 
surface  from  the  cylindrical  surface  is  relatively  small,  especially 
along  the  hub  and  casing  walls. 

The  following  steps  are  therefore  suggested: 

(a)  In  the  initial  calculation,  the  flow  surfaces  are  assumed  to  be 
cylindrical  and  the  set  of  equations  (60)  to  (65)  derived  for  cylindrical 
flow  or  the  approximate  method  given  in  reference  39  can  also  be  used  to 
obtain  the  streamlines  and  circumferential  variation  of  fluid  state  on 
Si  surfaces  at  three  or  more  .radii. 
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(Id)  From  the  data  obtained  in  step  (a),  an  S2  stream  surface  about 
midway  between  two  blades  is  constructed  by  connecting  the  streamlines 
which  divide  mass  flow  on  the  Sp  surfaces  in  the  same  percentages.  The 
direction  numbers  of  the  surface  and  the  Wu  and  Wz  at  the  surface 
obtained  in  step  (a)  give  the  starting  value  of  Wr  by  use  of  equa¬ 
tion  (31).  The  factor  B  is  evaluated  either  directly  from  the  angular 
distances  between  streamlines  obtained  in  (a)  or  according  to  equa¬ 
tion  ( 100a )  with  C  evaluated  from  the  information  obtained  in  (a).  Its 
value  at  other  radii  is  obtained  by  interpolation  or  by  proportioning 
according  to  the  channel-width  ratio.  Calculation  of  the  flow  on  this 
surface  is  then  made  by  the  use  of  equations  (9l)  to  (115).  For  subsonic 
flow  with  irrotational  inlet  flow,  the  solution  obtained  on  the 
surface  is  easily  extended  circumferentially  by  series  expansion  using 
equations  (117)  to  (130).  The  values  obtained  can  be  further  adjusted 
to  fit  the  given  blade  (reference  39)  and  can  be  used  in  a  more  accurate 
second  calculation  on  Sp  surfaces  in  the  next  step.  For  subsonic  flow 
with  large  rotationality  at  the  inlet  and  supersonic  flow  with  signifi¬ 
cant  check  caused  by  the  blade  entrance  angle,  it  is  more  desirable  to 
obtain  the  information  on  circumferential  variations  by  the  use  of  two 
or  more  S2  surfaces  at  or  near  the  two  blade  surfaces. 

(c)  The  radial  variation  of  fluid  state  computed  from  the  solution 
obtained  in  step  (b)  or  the  variation  of  the  radial  distance  between 
streamlines  is  used  to  determine  the  factor  b  and  used  in  the  principal 
equation  (48)  for  a  more  accurate  determination  of  Sp  surfaces  and  the 
flows  thereon.  The  general  equations  (32)  to-  (5l)  should  now  be  used 
for  the  Sp  surfaces  located  between  hub  and  casing,  if  not  at  or  near 
these  walls. 

(d)  The  calculation  of  S2  surfaces  can  again  be  repeated  and  so 
forth. 

If  the  inlet  flow  is  quite  rotational,  so  that  the  Sp  surfaces 
along  the  walls  and  the  S2  surface  near  the  blades  may  turn  around  the 
corners,  these  surfaces  should  be  chosen  at  a  short  distance  from  these 
boundary  walls  as  shown  in  figure  3.  By  the  use  of  these  two  kinds  of 
surface,  the  secondary  flow  caused  by  a  rotational  inlet  profile  or  by 
the  turning  of  the  blades  is  included  in  the  complete  solution. 

Axial  turbomachines  with  tapered  or  curved  walls.  -  The  steps 
involved  here  are  quite  similar  to  those  of  the  preceding  case,,  except 
that  for  the  initial  calculation  of  Sp  surfaces  along  or  near  the 
tapered  or  curved  wall,  either  equations  (52)  to  (59)  are  to  be  used,  or 
equations  (13* )  to  (23!)  given  in  reference  (35)  can  be  used  in  the 
manner  given  in  reference  (39). 
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Radial-  and  mixed-flow  type  turbomachines  with  curved  vails.  -  In 
this  type  of  machine  it  is  not  desirable  to  start  the  computation  on  the 
Sp  surfaces  because  the  flow  surfaces  near  the  walls  may  deviate  con¬ 
siderably  from  surfaces  of  revolution  because  of  the  long  flow  path.  On 
the  other  hand,  the  solidity  of  the  blade  is  very  high  and  the  blade 
section  is  uniformly  thin.  As  a  result,  the  shapes  of  the  Sg  surfaces 
are  closely  related  to  the  blade  shape  and  the  factor  B  can  be  esti¬ 
mated  relatively . accurately  from  the  blade  thickness  distribution.  The 
following  steps  are  therefore  suggested: 

(a)  The  computation  is  begun  on  the  'S2,m  surface.  For  subsonic 
irrotational  inlet  flow,  computation  need  be  made  only  on  a  mean  Sg 
surface  and  the  solution  can  be  extended  out  circumferentially  by  equa¬ 
tions  (117)  to  (130).  The  approximate  method  given  in  reference  40  can 
also  be  used  in  the  initial  calculation.  For  subsonic  flow  with  rota¬ 
tional  Inlet  profile  and  for  supersonic  flow  it  may  again  be  more 
desirable  to  compute  two  or  more  S2  surfaces  between  the  blades. 

(b)  The  data  obtained  in  step  (a)  may  be  used  to  make  calculations 
for  three  or  more  Sp  surfaces  between  hub  and  casing  walls. 

(c)  The  processes  (a)  and  (b)  can  be  repeated  until  the  desired 
accuracy  is  reached. 


Inverse  Problem 

Conditions  prescribed  on  mean  stream  surface.  -  In  the  inverse  or 
design  problem  it  is  most  convenient  to  consider  a  mean  stream  surface 
of  the  S2  kind  about  midway  between  two  neighboring  blades  to  be 
designed  (figs.  3  to  5).  From  the  results  developed  previously  for 
such  surfaces,  it  is  seen  that  in  addition  to  the  factor  B,  the 
designer  can  specify  only  one  relation  among  the  fluid  properties  on 
that  surface,  which  can  be  either  a  velocity  component,  a  relation 
between  two  velocity  components,  or  one  other  reasonable  condition. 

The  factor  B  essentially  controls  the  blade  thickness  distribution, 
whereas  the  relation  specified  on  the  surface  essentially  controls  the 
mean  camber  surface  of  the  blade.  From  a  consideration  of  strength  and 
Mach  number  in  general,  and  the  requirement  of  coolant  passage  in  the 
case  of  cooled  turbine  blades,  the  designer  always  has  a  very  good  idea 
of  what  kind  of  blade -thickness  distribution  he  wants.  With  this  thick¬ 
ness  distribution,  the  ratio  of  pitch  minus  circumferential  thickness 
of  blade  to  pitch  can  be  obtained.  After  correcting  this  ratio  with 
some  known  relations  between  this  ratio  and  B  ( such  as  those  given  in 
references  29,  30,  and  35),  especially  near  the  leading  and  trailing 
edges,  it  can  be  taken  as  the  factor  B  in  equation  (lOl).  Then  from 
the  type  of  velocity  diagram  or  a  certain  feature  of  blade  shape 
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desired,  a  relation  along  the  mean  stream  surface  Sg  m  can  be  pre¬ 
scribed  and  coordinates  of  the  mean  surface  and  the  flow  on  that  surface 
can  be  solved  at  the  same  time  by  equations  (lOl)  to  (115).  It  may  be 
noted  that  in  this  process,  the  designer  still  has,  in  general,  a  little 
freedom  in  choosing  the  value  of  z0  in  equation  (115c) .  For  a  rota¬ 
ting  blade,  z0  is  usually  taken  somewhere  near  the  center  of  gravity 
of  the  blade  section,  whereas  for  the  stationary  blade,  the  position  of 
z0  can  be  utilized  to  control  the  magnitude  and  distribution  of  Fr 
in  the  most  favorable  manner. 

Boundary  conditions  for  mean  stream  surface.  -  In  the  solution  of 
this  S2  m  surface,  the  boundary  conditions  are  a  little  different  for 
subsonic  and  supersonic  flow.  For  subsonic  flow,  not  only  the  varia¬ 
tions  of  the  stream  function  at  stations  far  upstream  and  downstream 
are  given,  the  meridional  contours  of  the  hub  and  casing  walls  are  also 
given  (these  contours  can  be  determined  by  approximate  calculations 
from  blade  row  to  blade  row  such  as  given  in  references  17  and  41) .  For 
supersonic  flow,  the  variation  of  the  stream  function  and  its  normal 
derivative  is  prescribed  on  an  initial  curve,  which  is  not  a  character¬ 
istic  curve.  Then  either  the  hub  and  casing  contours  are  prescribed, 
or  only  the  casing  contour  but  with  one  more  velocity  component  along 
the  casing  is  prescribed.  In  the  second  case,  the  flow  is  extended 
toward  the  axis  of  the  machine  and  the  hub  contour  is  determined  by  the 
shape  of  the  last  streamline  for  the  required  mass  flow. 

Determination  of  blade  shape.  -  For  subsonic  irrotational  flow, 
the  solution  obtained  on  the  mean  stream  surface  can  be  extended  out 
circumferentially  by  using  equations  (117)  to  (130).  The  blade  sur¬ 
face  can  be  then  determined  as  follows: 

(a)  The  position  of  the  mean  stream  surface  is  first  determined 
by  solving  the  circumferential  coordinate  as  a  function  of  the  axial 
coordinate  at  several  radii.  With  the  circumferential  coordinate 
measured  from  the  radial  element  of  the  surface  chosen  at  Zq ,  equa¬ 
tion  (116)  gives,  at  a  constant  r: 


z=zQ 

(b)  The  blade  coordinates  (r,  cp )  will  first  be  chosen  at  one 
station  z0  as  follows  (see  fig.  9):  The  mass  flow  passing  through 
the  z0  plane  between  the  mean  stream  surface  and  the  tentative  suction 
or  pressure  surface  is  computed  as  follows: 


Zo 


dx 


u/m 


(131) 
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frc 

p  Wzr  dr  dcp 

Jrh 


(132) 


J  rh 


p  Wz  r  dr  dCp 


(133) 


Because  of  the  inaccuracy  in  B  for  the  Hade -thickness  effect,  the 
mass  flow  obtained  will  he  a  little  different  from  that  required.  The 
blade  coordinates  Cf^  and  cpp  as  functions  of  zD  and  r  are  modi¬ 
fied  until  the  mass  flow  checks.  It  is  not  important  that  %  and  Mp 
are  a  little  different  from  one -half  the  required  mass  flow  as  long  as 
their  sum  is  equal  to  the  total  mass  flow,  but  once  the  division  is 
chosen,  it  should  be  maintained  at  other  z-stations. 

(c)  The  blade  coordinates  obtained  at  z  =  z0  are  extended 
upstream  and  downstream  according  to  the  velocity  components  evaluated 
at  the  surface.  For  example,  for  a  short  distance  z  -  zD  away,  .the 
changes  in  the  blade  surface  coordinates  r  and  cp  are 


(134) 


(135) 


After  r  and  Cp  are  thus  obtained,  the  total  mass  flow  may  be  checked 
again  by  equations  (l32)  and  (133). 

When  the  blade  coordinates  are  obtained  close  to  the  leading  and 
trailing  edges,  they  can  be  faired  in  according  to  some  standard  shapes. 
A  blade  shape  is  therefore  obtained  in  which  the  three-dimensional  flow 
of  the  fluid  is  considered.  The  right  kind  of  three-dimensional  blade- 
thickness  distribution  is  obtained  and  a  good  knowledge  about  the  flow 
on  the  blade  surface  is  also  ayailable  at  the  same  time.  The  data 
obtained  in  the  solution  can  also  be  used  directly  for  a  more  accurate 
and  detailed  determination  of  the  velocity  variation  around  the  nose  of 
the  blade,  for  a  relatively  quick  check  of  the  series  approximation,  or 
for  improvement,  if  necessary,  of  the  inverse  solution  throughout  by 
the  method  given  earlier  for  solving  the  direct  problem.  This  process 
seems  to  be  the  quickest  way  of  establishing  some  standard  three- 
dimensional  flow  variations  for  typical  designs  of  blades  from  which  a 


58 


MCA  TN  2604 


good  approximate  method  for  routine  design  calculations  can  he  estab¬ 
lished  and  of  providing  a  basis  on  which  the  viscous  flow  along  the 
blade  surfaces  and  hub  and  casing  walls  can  be  analyzed.  The  results 
given  in  references  29  and  35  indicate  that  for  blades  of  high 
solidity,  three  terms  in  the  series  give  sufficiently  accurate 
results. 

For  subsonic  flow  with  vorticity,  the  circumferential  extension 
cannot  be  accurately  made  at  the  present  because  of  lack  of  adequate 
methods  for  the  determination  of  vorticity  variation  along  the  mean 
stream  surface  An  estimate  of  this  variation  can  be  made, 

however,  and  the  solution  can  be  cheeked  later.  An  alternative  method 
is  to  use  the  shapes  of  the  streamlines  and  the  distances  between  them 
obtained  in  the  &2,m  solution  and  to  design  the  blades  with  the 
assumption  that  the  flow  surfaces  are  surfaces  of  revolution  by  the 
method  given  in  reference  35.  Inasmuch  as  the  rotationality  of  inlet 
flow  is  usually  serious  only  in  later  stages  of  a  multistage  compressor 
where  the  hub-tip  radius  ratio  is  high,  this  assumption  is  reasonably 
good. 


For  supersonic  flow,  the  flow  in  the  mean  stream  surface  Sg^m  is 
also  determined  first.  If  the  shock  due  to  the  entrance  wedge  angle 
is  small,  an  approximate  solution  of  the  blade  shape  can  also  be 
obtained  by  the  series  method  neglecting  the  finite  jump  across  the 
shock  or  using  an  estimated  value.  The  improvement  of  the  flow  varia¬ 
tion  for  the  resultant  blade  is  then  more  important  than  that  in  the 
subsonic  case.  Local  modification  of  the  blade  shape  can  also  be  made 
if  the  velocity  distribution  on  the  blade  obtained  is  unsatisfactory. 

A  better  method  is  to  use  the  shape  of  the  streamlines  and  the  dis¬ 
tances  between  them  obtained  in  the  Sg^m  solution  and  to  design  the 
blades  assuming  flow  surfaces  of  revolution  according  to  the  method 
given  in  reference  30. 

The  processes  described  here  for  the  three-dimensional  solution 
have  been  and  are  now  being  used  to  analyze  the  compressible  flow 
through  a  number  of  compressors  and  turbines.  Some  of  the  results 
obtained  are  given  in  reference  42. 


GENERAL  METHODS  OF  SOLVING  PRINCIPAL  EQUATION 

In  the  solution  of  the  Sp  surface  for  the  direct  problem  and  of 
the  Sg  surface  for  both  the  direct  and  inverse  problems >  the  main 
calculation  is  the  solution  of  the  principal  equation,  which  is  a 
second-order,  nonlinear  partial  differential  equation  in  two  independent 
variables.  The  case  when  the  principal  equation  is  elliptic  will  be 
considered  first. 
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Elliptic  Case 

A  common  form  of  the  principal  equation  is  written  as  follows: 

Mg (bp) 2  =  0  (136) 

In  equation  (136),  \|r  and  b  are  used  for  both  Sp  and  Sg  surfaces; 

T]  denotes  Cp  for  the  Sp  surface  and  r  for  the  S2  surface; 

£  denotes  z  or  r  for  the  Sp  surface,  and  z  for  the  Sg  sur¬ 
face.  The  values  of  rj,  t, ,  J,  and  K  for  each  individual  case  are  given 
in  the  following  table: 


iaj  1  t JL 
nan  ar2 


+  j 


ai 


+  K 


a2\|r 

an2 


+  L 


a£ 
an  ' 


Case 

Surface 

Coordinates 

j 

K 

L 

Equation 

n 

. 

1 

■ 

Sp  (general) 

cp 

z 

0 

Mm 

nu 

nrr 

(48a) 

p 

(surface  of 
revolution) 

cp 

K 

0 

(59a) 

Ci 

z 

r 

r2 

3 

Si  (cylindrical 
surface) 

cp 

•  z 

■■0 

! 

1 

r2 

0 

(65a) 

4 

Si  (general) 

Cp 

r 

0 

1 

r2 

0 

(75a) 

Si  (surface  of 
revolution) 

Cp 

1 

1  +x“2 

n 

(81a) 

•j 

r 

r 

r2 

6 

Si  (radial 
plane ) 

Cp 

r 

1 

r2 

0 

(86a) 

7 

s2  (vur  given) 

r 

z 

0 

1 

_  1 

r 

(107a) 

8 

s2  specified^ 

r 

z 

0 

1  +  g2 

_  i  +  g 

r  &  3r 

(114a) 
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The  equation  is  nonlinear  even  in  the  case  of  incompressible  flow.  In 
the  numerical  computation,  it  is  convenient  to  rewrite  the  equation  in 
the  following  form: 


where 


+  K 


a2\|r 


a  In  bp  3jr  a  In  bp  gjf 

ari  at]  at;  at; 


Mxbp 


M2(~bp)2 

aj 

3T) 


(137) 


(138) 


and  is  evaluated  from  any  approximate  solution  at  the  start  of  the  cal¬ 
culation  and  from  the  values  of  i|r  and  p  obtained  in  the  previous 
cycle  during  the  calculation.  For  simple  boundary  shapes  for  an  S2 
surface  and  simple  functions  of  J,  K,  and  L,  it  is  possible  to  find 
a  Green's  function  G(t),.£,  X.,  y)  with  its  proper  characteristics  so 
that  the  solution  of  the  problem  can  be  written  in  the  following  form 
(for  example,  see  reference  10): 


t  (IjO  =  SS  G  C1);  xj  y)  N  (x>  y)  clx  dy  (139) 

If  the  boundary  wall  is  arbitrarily  curved,  it  is  necessary  in  this 
method  to  use  the  technique  of  conformal  transformation  to  render  the 
given  boundary  into  a  simpler  one,  such  as  cylindrical.  Because  this 
process  will  involve  a  numerical  solution  of  the  Laplace  equation  with 
the  given  boundary  shape,  it  may  be  better  to  solve  the  given  equa¬ 
tion  (137)  directly  with  the  given  shape  by  the  numerical  method. 
Furthermore,  this  method  will  be  the  only  choice  in  the  general  case 
where  J  and  K  are  complicated  functions,  which  makes  the  task  of 
obtaining  the  proper  Green1 s  function  very  difficult  if  not  impossible. 

Finite-difference  form  of  principal  equation.  -  In  order  to  solve 
the  given  equation  (137)  directly,  the  general  numerical  differentiation 
formula  for  first  and  second  derivatives  with  the  function  value  given  at 
unequally  spaced  grid  points  using  second-  and  higher-degree  polynomial 
representation  as  given  by  reference  26  is  used  to  give  the  finite- 
difference  expressions  conveniently  and  accurately  at  the  grid  point 
near  the  curved  boundary.  If  the  value  of  any  quantity  q  on  the 
stream  surface  under  investigation  corresponding  to  a  number  of  values 
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of  one  of  the  independent  variables  x  not  equally  spaced,  denoted  by 

x°,  x1,  .  .  .  x1,  .  .  .  xn,  is  given,  the  m^*1  derivative  of  q  (on 
the  surface)  with  respect  to  x  when  x  =  x*  may  be  written 


O^q) 


x=x1 


y  V  +  v 

i—t  n  j  n 

j=0 


(140) 


The  differentiation  coefficients  B  and  the  coefficients  of  the  deriva¬ 
tives  in  the  first  or  second  remainder  term  have  been  explicitly 
expressed  in  reference  26  in  terms  of  the  spacings  between  the  successive 
grid  points  by  using  polynomials  of  the  second,  third,  and  fourth  degree 
for  general  nonuniform  spacing  throughout  and  for  the  special  case  near 
a  tapered  or  curved  boundary  where  only  the  first  or  last  spacing  is 
different  from  the  others.  For  the  special  case,  these  coefficients 
have  also  been  computed  for  different  ratios  of  the  distance  between  the 
boundary  and  the  nearest  point  and  the  other  spacing,  from  0.1  to  1.29 
in  intervals  of  0.01,  and  are  given  in  reference  26.  For  spacing  ratios 
lying  between  these  tabulated  intervals,  B  can  be  obtained  from  the 
values  tabulated  by  applying  interpolation  formulas  given  in  refer¬ 
ence  43,  or  by  the  direct  use  of  the  formulas.  Differentiation  coeffi¬ 
cients  B  for  equal  intervals  using  various  degrees  of  polynomials  are 
given  earlier  by  Bickley  in  reference  44. 


In  the  present  fluid-flow  problems,  a  large  region  must  be  covered 
in  order  to  get  to  the  boundary  conditions  which  are  always  given  at 
stations  far  upstream  and  downstream  of  the  blade  row.  In  order  to 
reduce  the  labor  of  computation,  it  is  desirable  to.  attempt  to  reduce 
the  number  of  grid  points  required  for  a  given  accuracy  by  using  a 
degree  of  polynomial  higher  than  the  customary  second.  A  study  of  the 
expression  of  the  remainder  terms  (see  reference  26)  and  actual  experi¬ 
ence  in  the  present  problem  show  that,  in  most  cases,  the  use  of  the 
fourth-order  polynomial  will  reduce  the  necessary  number  of  grid  points 
to  less  than  one-fourth  that  required  by  the  second-order  polynomial. 
Near  the  leading  and  trailing  edges  of  the  flow  on  surfaces  of  the  Sq 
kind,  the  variation  of  i|r  is  such  that  accuracy  is  obtained  most 
effectively  by  using  small  spacing  there.  In  such  case,  the  grid 
pattern  should  be  chosen  at  these  regions  first,  and  either  be  kept 
constant  or  be  continuously  increased  toward  the  inlet  and  exit  stations. 


With  the  grid  pattern  and  the  order  of  polynomial  representation 
selected,  the  coefficients  B  at  each  point  can  be  obtained  from  refer¬ 
ences  44  and  26  for  equally  and  unequally  spaced  points.  Then  the  dif¬ 
ferential  equation  (137)  at  any  grid  point  whose  i|r  value  is  is 

replaced  by  the  following  algebraic  equation: 
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* 

(2A*  +  xAi)  i(rk  -  N1  =  0  (141) 

vn  j  n  j'  Y  v 

on  the  surface  considered 

corresponding  to  the  grid  points  along  constant  £  and  constant  t) 
lines,  respectively.  (See  figs.  10  and  11.)  It  should  he  noted  that, 
in  accordance  with  the  definition  of  the  special  partial  derivatives, 

\|r  values  are  those  on  the  surface  S;  whereas  the'  grid  spacings 
involved  are  just  the  distances  along  the  tj-  and  ^-coordinates. 

Boundary  conditions.  -  In  flow  on  surfaces  of  the  first  kind,  the 
flow  picture  is  as  shown  in  figure  9.  Arbitrarily  assigning  a  value 
\|fj  on  the  suction  surface,  the  value  on  the  pressure  surface  of 

the  next  blade  is  determined  from  the  mass  flow  passing  between  them. 

These  two  values  are  used  as  the  end  values  in  equation  (l4l)  for  grid 
points  next  to  the  boundary.  Outside  the  blade  region,  however,  the 
position  of  the  dividing  streamline  is  not  known.  Instead,  there  is 
the  condition  that  the  flow  repeats  itself  or  the  i|r-value  increases 
■by  ^n-^i  when  cp  increases  by  an  amount  equal  to  the  pitch  angle 

(2it  divided  by  number  of  blades).  It  is  then  convenient  to  draw  any 
two  parallel  lines  up  to  the  leading  and  trailing  edges  of  the  blade  and 
consider  only  the  grid  points  lying  between  the  two  reference  lines. 

For  the  cp- derivative  at  a  point  tyc,  for  example,  the  required  i|rb 
value  is  obtained  from  >|rf,  which  is  a  pitch  angle  away  from 
(fig.  10),  as 


i  (J1  nAi  +  K±  nAi)  +  S 
J=0  J  J  k=0 


where  ijr1  and  ^  denote  the  values  of 


=  *f  -  (tIZ  -  *j)  (142) 

This  relation  is  used  between  the  inlet  station  1-1  and  the  leading  edge 
of  the  blade  and  between  the  exit  station  2-2  and  the  trailing  edge  of 
the  blade  when  the  Sq  surfaces  are  assumed  to  be  surfaces  of  revolu¬ 
tion.  For  the  general  Sq  surface  where  its  deviation  from  the  surface 
of  revolution  is  considered,  modification  has  to  be  made  in  places  such 
as  shown  in  the  exit  portion  of  figure  9.  Because  of  the  twist  of  the 
flow  surface,  the  dividing  line  from  station  1-1  to  the  leading  edge  of 
the  blade  becomes  two  separate  lines  from  the  trailing  edge  of  the  blade 
to  the  exit  station  2-2,  accompanied  by  trailing  vortices.  Although  the 
flow  still  repeats  itself  circumferentially  every  pitch  angle,  the  use 
of  equation  (142)  for  the  derivative  at  a  point  close  to  these  lines 
will  give  inaccurate  results.  In  these  cases,  it  is  better  to  use  the 
end-point  differentiation  formulas  to  evaluate  the  derivatives. 
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At  the  station  1-1  sufficiently  far  upstream  of  the  blades,  the 
flow  condition  can  be  taken  as  uniform  and  the  flow  angle,  equal  to  the 
given  inlet  angle.  For  the  point  h,  the  \|r  value  at  point  i 
upstream  can  be  obtained  from  the  given  flow  angle  as  follows: 


t1  = 


5^  tan 


-  *€) 


=  (l  -  jp-  tan  \|r®  +  tan  llr*1 


(143) 


Thus,  the  required  value  upstream  of  station  1-1  can  be  replaced  by 
the  values  on  that  station,  and  only  the  \[r  values  downstream  of 
station  1-1  will  be  involved  in  the 'finite-difference  expression  (l4l). 


An  alternative  method  to  take  account  of  the  inlet  condition  is 
as  follows:  If  the  first  station  1-1  is  chosen  sufficiently  far  from 
the  blades,  the  variation  of  the  stream  function  upstream  of  the 
station  1-1  is  linear  in  the  circumferential  direction.  The  value  of 
the  stream  function,  however,  depends  on  the  inlet  angle.  If  solutions 
for  a  range  of  inlet  angle  are  desired,  they  can  be  obtained  by  speci¬ 
fying  a  number  of  sets  of  linearly  varying  stream  functions  upstream  of 
station  1-1  as  fixed  boundary  values.  The  slope  of  the  streamlines 
obtained  in  the  solution  at  the  inlet  then  gives  the  value  of  the  inlet 
angle.  If,  however,  the  solution  for  a  certain  specific  inlet  angle 
is  desired,  the  streamline  obtained  in  the  solution  must  be  adjusted 
according  to  that  inlet  angle,  for  example,  as  jk  in  figure  10  is 
adjusted  to  position  .  gk,  thereby  obtaining  an  improved  set  of  boundary 
values  of  the  stream  functions  to  be  used  in  the  next  calculation. 

This  method  is,  of  course,  not  so  accurate  and  convenient  as  the  previ¬ 
ous  method  for  obtaining  a  solution  for  a  given  inlet  angle,  but  is 
desirable  in  the  matrix  solution  because  the  inlet  angle  is  then  not 
involved  in  the  matrix  factorization,  thereby  making  the  same  matrix 
factors  usable  for  a  range  of  inlet  angle  and  Mach  number. 


At  the  exit  station  far  downstream'  of  the  blade,  the  same  methods 
can  be  applied.  For  a  blade  having  a  sharp  trailing  edge,  the  Kutta- 
Joukowski  condition  can  be  used  and  the  correct  exit  angle  far  down¬ 
stream  is  the  one  that  gives  the  flow  at  the  trailing  edge  satisfying 
that  condition.  For  round  trailing  edges,  either  the  position  of  the 
stagnation  point  is  assumed  or  some  available  empirical  rule  for  the 
exit  angle  is  used.  If  the  calculation  is  made  to  compare  with  certain 
experimental  results,  the  measured  exit  angle  may  be  used. 
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In  flow  along  surfaces  of  the  second  kind,  the  boundary  walls 
extend  all  the  way  to  the  inlet  and  exit  stations  with  the  \|f  values 
given  on  the  walls  (fig.  11 ) .  Across  the  inlet  and  exit  stations,  the 
flow  is  considered  to  he  uniform  and  parallel  to  the  walls  so  that  the 
required  i|r  value  outside  the  station  can  he  obtained  by  an  equation 
similar  to  equation  (143) .  For  the  inlet  station  where  the  axial 
velocity  is  radially  uniform  and  there  is  no  radial  or  tangential 
velocity,  \|r  varies  as  the  square  of  the  radius.  For  the  exit  station 
with  a  certain  radial  gradient  in  fluid  state,  the  radial  variation  of 
♦  can  be  determined  from  the  corresponding  radial  variation  in  axial 
velocity  and  density. 

Solution  of  finite-difference  equations.  -  With  the  grid  system 
and  the  degree  of  polynomial  representation  chosen  and  the  boundary 
conditions  taken  into  account,  the  problem  remaining  is  the  solution 
of  the  set  of  linear  algebraic  equations  (l4l)  written  for  all  interior 
grid  points.  For  a  small  number  of  solutions  with  a  given  blade,  the 
best  method  is  the  relaxation  method  (references  25,  33,  45,  and  36). 

A  modification  of  this  method  involving  the  use  of  higher-order  differ¬ 
ences  is  suggested  by  Fox  (reference  46).  Formulas  and  tables  of 
coefficients  obtained  in  reference  26  enable  the  direct  use  of  higher- 
degree  polynomials  for  problems  with  curved  boundaries  (reference  29). 
For  the  present  flow  problems,  it  is  necessary  to  include  a  large 
domain  to  get  to  the  boundary  conditions  that  are  given  at  places  far 
from  the  blades,  and  the  use  of  higher-degree  polynomials  whenever  it 
is  applicable  greatly  reduces  the  numerical  work. 

If  a  number  of  cases  are  to  be  solved  for  a  given  geometry  (same 
blades  for  S±  surface  and  same  hub  and  casing  shapes  for  S2  sur¬ 
face),  it  is  advantageous  to  solve  the  problem  on  a  large-scale  digital 
computing  machine.  If  a  high-speed  digital  machine  is  available,  the 
simultaneous  equations  may  be  solved  by  Liebmannfs  iterative  process, 
which  is  the  most  simple  to  set  up.  For  quicker  results  or  when  only  a 
relatively  slow- speed  machine  is  available,  the  matrix  process  discussed 
in  reference  26  is  most  suitable.  In  a  calculation  of  the  S2,m  sur¬ 
face  for  a  gas  turbine  and  in  a  calculation  of  the  S±  surface  of 
revolution  for  a  centrifugal  compressor,  the  coefficient  matrices 
(about  400  and  200  interior  grid  points  for  the  two  problems,  respec¬ 
tively,  and  the  fourth-degree  differentiation  formula  are  used)  were 
factorized  into  the  lower  and  upper  triangular  matrices  on  an  IBM  CFEC 
and  an  IBM  604,  respectively,  in  about  60  hours.  The  determination  of 
t  for  a  given  set  of  values  of  N  took  only  2  hours  on  the  CFEC  for 
the  gas-turbine  problem.  The  gas-turbine  problem  was  also  worked  out 
on  an  Univacj  the  factorization  took  only  11  minutes  and  the  determina¬ 
tion  of  )|r,  2.5  minutes.  The  increasing  availability  of  these  high¬ 
speed  large-scale  digital  calculating  machines  will  render  the  suggested 
method  of  solving  the  three-dimensional-flow  problem  a  practical  one. 
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General  table  for  evaluation  of  density  from  ^-derivatives.  - 
After  the  \[f  values  are  obtained  at  the  end  of  each  cycle  of  calcula¬ 
tion,  the  velocity  components  are  evaluated  from  the  derivatives  of 
with  respect  to  the  coordinates,  after  the  density  is  obtained  as 
follows:  From  equations  (46),  (57),  (74),  (80),  (85),  (lG5),  and  (113), 
the  relation  between  h  or  p  and  ^-derivatives  can  be  put  into  a 
common-  form  as 

h  ,  j  +  <°|!  .  X  .  i  (*p)-2  [tt  (J  -ft)2  +  (i  ft)2]  (144) 


The  quantities  represented  by  X,  as  well  as  by  €,  t],  and  £  for  dif¬ 
ferent  cases,  are  given  in  the  following  table: 


With  the  ^-derivatives  evaluated,  if  an  exact  determination  of  h 
or  p  from  the  preceding  equation  considering  the  variation  of  specific 
heat  with  temperature  is  desired,  the  Keenan  and  Kay  gas  tables  (refer¬ 
ence  47)  can  be  used.  With  two  or  three  readings  of  h  and  p  (or  its 
reciprocal,  specific  volume),  the  correct  value  of  h  or  p  satisfying 
equation  (144)  is  found.  For  most  cases  where  the  temperature  range 
involved  is  not  too  large  the  use  of  an  appropriate  average  value  of  y, 
y,  may  give  accurate  enough  results.  With  the  use  of  an  average  y, 
the  density  at  any  point  in  the  flow  field  can  be  related  to  the  inlet 
total  value  by  equation  (12a)  as 
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T-l  s$  ±-s*  I  + 


a/'r2 


-*  [4sT*(Jstl 


1 

2lV-l 


2("bp)  % 


*  * 

sTji-s 


(145) 


In  order  to  make  out  a  general  table  for  the  calculation  of  density  from 
the  ^-derivatives,  the  preceding  equation  is  rewritten  as 

2 


■  (■  -  s 


(146) 


where 


I  +  i  co2r2  -  X  \  ^  s*-s* 


VO- 

*  -  [K  (f  If)2  +  I  f)2](2Ei)  1  <b»T,i> 


-2/l  +  |o>2r2  -  A1"'1  2(s*-s*  ) 

V  % — 7  e 


The  functional  relations  between  Z  and  $  are  given  in  table  I  for  y 
equal  to  1.4  and  4/3,  respectively.  From  the  given  inlet  condition  and 
the  given  X  values,  the  variation  of 

/  2  2  \  T-l 


(2  Hi)”1(hp,p^i) 


is  first  competed  and  plotted  as  an  auxiliary  graph  or  table  as  a  func 

tion  of  - X  ^  .  A  similar  auxiliary  graph  or  table  is  prepared 

for  the  variation  of 


122  \  T-l 

I  +  |co2r2  -X] 


as  a  function  of 

[  (S?X^ 

value  of  1  -~— 


/o)2r2  \ 

^ _ —  -  X  J.  Anytime  during  the  calculation,  from  the 
-  X  j  at  each  point  (in  general,  X  changes  during 


successive  improvements  between  S-]_  and  Sg  surfaces), 
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(2Hi)"1(bpT>1)-2 


I  + 


\  o2r2  -  x\ 

~T—) 


_  T+1 

r-i 


is  read  from  the  first  graph  or  table  and  is  combined  with 

k(  I  §tf  +fi  ®ix2 

k  \r  at]/  +\e  9^ 

and  the  entropy  factor  to  obtain  $  .  The  value  of  2  is  then  read  from 
tables  I  or  II.  After  the  value  of 


I  +  |  o>2r2  -  X 

% 


1 

r-i 


is  read  from  the  second  curve  or  table,  the  density  ratio  is  obtained. 


Hyperbolic  Case 


In  the  hyperbolic  case,  the  main  problem  is  the  solution  of  the 
following  principal  equation,  written  in  a  common  form  for  the  two  kinds 
of  flow  surface: 


3^ 


2  £  2?±_ 
V  9^3r) 


+  4^|  +  m 


at  ,  n  at 

3^  +  V  3T] 


=  o 


(147) 


with  the  initial  condition  that  t  and  its  normal  derivative  are  given 
on  a  curve  which  is  not  a  characteristic  curve.  From  equation  (147) , 
the  equation  of  the  characteristic  curve  is 


J 


The  slopes  of  the  characteristic  curves  are 


(148) 


Al~  ■  (149a) 
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Lz-  +  (149b> 

The  coefficients  J,  K,  L,  and  v  and  the  independent  variables 
r)  and  £  for  the  eight  cases  considered  are  given  in  the  table  on  the 
following  page.  Using  these  values  of  J ,  K ,  and  L,  Ap  and  Ag  are 

also  expressed  in  terms  of  the  velocity  components.  Except  for  cases  2, 
5,  and  8,  they  can  also  be  expressed  in  the  usual  trigonometric  form, 
tan  (X  i  M-) .  The  values  of  X  and  \x  are  also  given. 

Changes  of  ^-derivatives  along  characteristic  curve.  -  When  the 
reference  point  on  the  -plane  moves  &long  the  image  of  the  charac¬ 
teristic  curve  in  the  rj^-plane  corresponding  to  a  small  change  in 

d£,  the  change  in  t)  is  drj  =  ^  d£.  Because  of  these  two  small 

changes,  the  change  of  any  quantity  q  on  the  surface  is  (fig.  12) 


dq 


d£ 


3t) 


da  _  §a  ,  A  aq 
d£  v  arj 


(150) 

(151) 


Hence  along 


d  aty  _  _a_  ajf  Ai  a  3i|r  32t|t  ^i  a2\|r 

df  at  ~  at;  a£  v  3ti  at  a^2  v 


(152) 


d  ay  _  3  3ty  +  ^-1  3  3\)f  _  a2i|f  +  ^1  a2\|f 

d£  9T]  at;  3t]  v  3^  3t|  9^3^  V  2 


(153) 


From  equations  (152)  and  (153), 


a2i|f  =  A  §!i 

a^n  dt;  an  “  v  8Tj2 


(154) 


_  A  §i  _d_  / ^i )  a^jf 

3^2  d£  at;  u  dt;  \v/  0^2 


(155) 
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Substituting  equations  (154)  and  (155)  into  equation  (147)  gives 


J  ^  If  +  (2K  -  jAi)  £  ^  f|  +  (JA1 


1 

2KAi  +  L)  -s-  — 
V2  3t] 


m  8ltf  4.  N  3l|r  n 

2+M9^+V8n=0 


(156) 


By  virtue  of  equations  (l49a)  and  (l49b),  equation  (156)  becomes 


d£  85  v  3t)  j  9^  j»  at) 

§ 

Similarly,  along  the  second  characteristic  curve  Ag, 


(157a) 


d  3f  4  d  3\|r  M  3*  IT  3t|t  _  - 

d£  95  V  d5  3r|  J  35  Jv  5i) 


(157b) 


Starting  from  two  points  a  and  b  a  short  distance  apart  on  the 
initial  curve,  equations  (149a)  and  (l49b)  give  the  tangent  to  the 
characteristic  curves  at  these  two  points  and  equations  (157a)  and 
(157b)  give  the  new  value  of  ch|r/5£  and  cty/dT]  at  the  point  of  inter¬ 
section  C  of  the  two  tangent  lines  (fig.  11 ) .  The  auxiliary  equa¬ 
tions  corresponding  to  the  particular  problem  are  then  used  to  deter¬ 
mine  other  pertinent  quantities  at  the  point  C.  This  process  is  to 
be  carried  step-by-step  downstream. 

Changes  of  fluid  velocity  and  direction  along  characteristic 
curve.  -  When  the  characteristic  curve  hits  the  boundary  wall,  it  is 
more  convenient  to  express  equations  (157a)  and  (157b)  in  terms  of  the 
magnitude  of  the  fluid  velocity  and  the  flow  direction.  In  order  to 
do  this,  the  definitions  of  ^-derivatives  are  first  put  in  a  common 
form  for  all  cases  as 

1^-  =  bp  r  (158) 

||  =  -  bp  €  (159) 

where  €  equals  1  and  r  for  the  and  S->  surfaces,  respectively. 

By  the  use  of  equation  (45), 
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d  Q\ 

d£  81 


TT  ,  i  d  In  b  .  1  dh 

--bpW,€  ^_3r-+_a?.. 


ds  *  d  l*1  "Wrj  d  In  £ 


(160) 


d  8*  TT  l  d  In  b  .  1  dh  ds *  .  d  ln  wt  .  d  In  r\ 


(161) 


Substituting  equations  (160)  and  (l6l)  into  equation  (156)  yields 


(a2  w5  -  w„)  (i-g 


d  ^  +  jl _  ^!V+  L  _  Ea\  + 

..  d£  a2d5  d£  )  +  V 2  a-z) 


A  TJ  d  111  r  TT  d  In  c  .  tt  N  rr  M  _ 

a2  w£  az  ~  wn  dC  w£  j  "  Wti  j  ~  0 


(162) 


¥^  =  w  sin  X 
Wf-  =  v  cos  X 


(163) 


cs2r2  1  /TI  2 


h  =  I  +  =~-  -  f  (W^  +  ¥^  +  ¥^)  = 


o>2r2 


|  (w2  +  Wfi2) 


(164) 


where  ¥|  is  equal  to  ¥r,  ¥r,  0,  ¥z,  ¥z,  0,  ¥u,  and  ¥u  for  cases  1 

to  8,  respectively.  By  the  use  of  equations  (163),  (164),  and  (144), 
equation  (l62)  can  be  written 


/  w£  \  dw 
\a2  "  /  d^ 


dw  +  cos  X  +  A2  sin  X  dx  _1_  _d_  /  +  <p2r2  wfc2\ 


Ao  cos  x  -  sin  x  d£  a2  d£ 


2  2 


d  In  Id 


2  cos  x  - 


d  In  6 


d  In  r\  ds 


sin  X  \finX  dg  "  ^  COsX~dT" J  +  d T  + 


1  M  sin  X  -  N  cos  X 
J  Ap  cos  X  -  sin  X 


(165) 
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A  similar  expression  can  be  obtained  for  the  change  in  w  and  X 
along  A2  by  replacing  A2  by  Aj_  in  the  preceding  equation.  For 
cases  1,  3,  4,  6,  and  1 ,  A  can  be  written  as  tan  (x  £  |i),  where  p.  is 

equal  to  sin-J-  through  purely  trigonometric  transformations  as 
w 

follows : 


A  K  +  /^K^-JL  ^  ±  a  /^w2-a2 

=  7  =  ’  a2  -  W^2 


sin  X  cos  X  ±  sin  p  cosp  _  sin  2x  ±  sin  2p 
2  cos  2X  +  cos  2{i 


cos^  X 


sm 


_  sin  (Xjhji)  cos  (X - )  cos  (x+p)  sin  (x-p) 
“  cos  (X+p)  cos  (X-p)  0r  COS  (X+p)  cos  (x-p) 


=  tan  (X+p)  or  tan  (X-p) 

For  these  cases,  equation  (165)  can  then  "be  written  (compare  refer¬ 
ence  30): 


1 

w 


dw  —  , 

-t-t  +  tan 

<3-^ 


dX 

c£ 


+  tan^ 


d  In  b  ,  ds*  ,  1  /  .  „  d  In  t 

d£  d£  A  cos  X  -  sin  X  lSin  X  d£ 


A  cos 


X 


d  In  r\ 
^  ) 


1  M  sin  X  -  N  cos  X  _  n 

J  A  cosX  -  sin  X  J 


(166) 


where  the  minus  and  plus  signs  on  the  second  term  and  subscripts  2 
and  1  for  X  in  the  last  two  terms  are  used  along  characteristics  Ai 

and  Ag,  respectively.  Equations  (165)  and  (166)  are  most  useful  when 

the  characteristic  hits  the  boundary  wall.  For  a  direct  problem,  the 
slope  there  is  known  from  the  given  blade  shape  and  for  an  inverse  or 
design  problem,  either  the  desired  turning  at  the  boundary  or  the 
velocity  on  the  boundary  is  prescribed.  With  either  dX  or  dw  known, 
dw  or  dX  is  evaluated  from  equation  (162)  or  (165)  (only  one  charac¬ 
teristic  equation  is  used  at  the  wall) .  For  convenience  of  setup  in 


v 
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calculation,  this  system  can  also  he  used  for  interior  points.  ■  Except 
that  more  terms  are  involved  in  the  -present  problem  and  that  w  takes 
different  meanings  in  different  cases,  the  procedure  of  calculation  is 
very  much  the  same  as  ordinary  two-dimensional  flow  described  in  refers 
ences  32  and  30. 


CONCLUDING  REMARKS 

A  general  theory  of  steady  three-dimensional  flow  of  a  nonviscous 
fluid  in  subsonic  and  supersonic  turbomachines  having  arbitrary  hub  and 
casing  shapes  and  a  finite  number  of  thick  blades  is  presented.  The 
solution  of  the  three-dimensional  direct  and  inverse  problem  is 
obtained  by  investigating  a  combination  of  flows  on  relative  stream 
surfaces  whose  intersection  with  a  z -plane  either  upstream  of  or  some¬ 
where  inside  the  blade  row  form  a  circular  arc  or  a  radial  line.  The 
equations  obtained  to  describe  the  fluid  flow  on  these  stream  surfaces 
show  clearly  the  several  approximations  involved  in  ordinary  two- 
dimensional  treatments.  They  also  lead  to  a  solution  of  the  three- 
dimensional  problem  in  a  mathematically  two-dimensional  manner  through 
an  iterative  process.  The  equation  of  continuity  is  combined  with  the 
equation  of  motion  in  either  the  tangential  or  the  radial  direction 
through  the  use  of  a  stream  function  defined  on  the  surface,  and  the 
resulting  equation  is  chosen  as  the  principal  equation  for  such  flows. 
The  character  of  this  equation  depends  on  the  relative  magnitude  of 
the  local  velocity  of  sound  and  a  certain  combination  of  velocity  com¬ 
ponents  of  the  fluid.  A  general  method  to  solve  this  equation  by  both 
hand  and  machine  computations  when  the  equation  is  elliptic  or  hyper¬ 
bolic  is  described.  The  theory  is  applicable  to  both  irrotational  and 
rotational  absolute  flow  at  the  inlet  of  the  blade  row  and  to  both 
design  and  off -design  operations. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  July  13,  1951 
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TABLE  I  -  GENERAL  DENSITY  TABLE 
(a)  f  =  1.4 


2 


5qp 

1 

Z 

2 

0.001 

.002 

.003 

.004 

.005 

1.0005007 

1.0010027 

1.0015061 

1.0020109 

1.0025171 

0.99949955 

.99899830 

.99849616 

.99799314 

.99748922 

0.006 

.007 

.008 

.009 

.010 

1.0030246 

1.0035335 

1.0040438 

1.0045555 

1.0050687 

0.99698452 

.99647894 

.99597249 

.99546516 

.99495686 

0.011 

.012 

.013 

.014 

.015 

1.0055833 

1.0060993 

1.0066168 

1.0071357 

1.0076561 

0.99444770 

.99393768 

.99342669 

.99291486 

.99240207 

0.016 

.017 

.018 

.019 

.020 

1.0081779 

1.0087012 

1.0092260 

1.0097522 

1.0102800 

0.99188844 

.99137386 

.99085834 

.99034199 

.98982460 

0.021 

.022 

.023 

.024 

.025 

1.0108093 

1.0113401 

1.0118724 

1.0124063 

1.0129417 

0.98930629 

.98878706 

.98826690 

.98774573 

.98722365 

0.026 

.027 

.028 

.029 

.030 

1.0134786 

1.0140171 

1.0145572 

1.0150989 

1.0156421 

0.98670066 

.98617666 

.98565167 

.98512569 

.98459881 

0.031 

.032 

.033 

.034 

.035 

1.0161869 

1.0167333 

1.0172814 

1.0178311 

1.0183825 

0.98407094 

.98354210 

.98301217 

.98248128 

.98194932 

0.036 

.037 

.038 

.039 

.040 

1.0189354 

1.0194900 

1.0200462 

1.0206042 

1.0211638 

0.98141649 

.98088260 

.98034775 

.97981176 

.97927482 

0.041 

.042 

.043 

.044 

.045 

1.0217252 

1.0222882 

1.0228530 

1.0234194 

1.0239876 

0.97873675 

.97819773 

.97765759 

.97711652 

.97657433 

0.046 

.047 

.048 

.049 

.050 

1.0245575 

1.0251292 

1.0257027 

1.0262780 

1.0268550 

0.97603112 

.97548680 

.97494137 

.97439485 

.97384733 

0.051 

.052 

.053 

.054 

.055 

1.0274338 

1.0280145 

1.0285970 

1.0291813 

1.0297674 

0.97329872 

.97274893 

.97219805 

.97164610 

.97109308 

5cp 

1 

Z 

Z 

0.056 

.057 

.058 

.059 

.060 

1.0303554 

1.0309453 

1.0315371 

1.0321307 

1.0327263 

0.97053890 

.96998357 

.96942708 

.96886954 

.96831077 

0.061 

.062 

.063 

.064 

.065 

1.0333238 

1.0339232 

1.0345246 

1.0351279 

1.0357332 

0.96775086 

.96718983 

.96662757 

.96606419 

.96549961 

0.066 

.067 

.068 

.069 

.070 

1.0363404 

1.0369496 

1.0375609 

1.0381742 

1.0387895 

.96493392 

.96436702 

.96379885 

.96322948 

.96265894 

0.071 

.072 

.073 

.074 

.075 

1.0394069 

1.0400263 

1.0406478 

1.0412714 

1.0418971 

0.96208713 

.96151415 

.96093991 

.96036442 

.95978768 

0.076 

.077 

.078 

.079 

.080 

1.0425249 

1.0431548 

1.0437869 

1.0444212 

1.0450576 

0.95920970 

.95863049 

.95804996 

.95746812 

.95688506 

0.081 

.082 

.083 

.084 

.085 

1.0456962 

1.0463370 

1.0469800 

1.0476253 

1.0482729 

0.95630069 

.95571503 

.95512808 

.95453976 

.95395006 

0.086 

.087 

.088 

.089 

.090 

1.0489227 

1.0495748 

1.0502292 

1.0508859 

1.0515449 

0.95335910 

.95276678 

.95217311 

.95157809 

.95098174 

0.091 

.092 

.093 

.094 

.095 

1.0522063 

1.0528701 

1.0535362 

1.0542048 

1.0548758 

0.95038397 

.94978478 

.94918428 

.94858229 

.94797890 

0.096 

.097 

.098 

.099 

.100 

1.0555492 

1.0562251 

1.0569034 

1.0575842 

1.0582675 

0.94737413 

.94676788 

.94616026 

.94555119 

.94494067 

0.101 

.102 

.103 

.104 

.105 

1.0589534 

1.0596418 

1.0603327 

1.0610263 

1.0617225 

0.94432862 

.94371513 

.94310022 

.94248371 

.94186569 

0.106 

.107 

.108 

.109 

.110 

1.0624213 

1.0631227 

1.0638268 

1.0645336 

1.0652431 

0.94124619 

.94062520 

.94000264 

.93937852 

.93875285 

5cp 

1 

Z 

Z 

0.111 

.112 

.113 

.114 

.115 

1.0659554 

1.0666704 

1 .0673882 
1.0681087 
1.0688320 

0.93812555 

.93749672 

.93686627 

.93623430 

.93560073 

0.116 

.117 

.118 

.119 

.120 

1.0695582 

1.0702883 

1.0710192 

1.0717540 

1.0724918 

0.93496548 

.93432857 

.93369008 

.93304993 

.93240806 

0.121 

.122 

.123 

.124 

.125 

1.0732326 

1.0739763 

1.0747230 

1.0754728 

1.0762256 

0.93176447 

.93111924 

.93047232 

.92982361 

.92917321 

0.126 

.127 

.128 

.129 

.130 

1.0769815 

1.0777405 

1.0785026 

1.0792678 

1.0800362 

0.92852106 

.92786714 

.92721149 

.92655410 

.92589489 

0.131 

.132 

.133 

.134 

.135 

1.0808079 

1,0815828 

1.0823610 

1.0831424 

1.0839272 

0.92523380 

.92457092 

.92390616 

.92323964 

.92257118 

.136 

.137 

.138 

.139 

.140 

1,0847153 

1.0855068 

1.0863017 

1.0871000 

1.0879018 

.92190089 

.92122868 

.92055458 

.91987858 

.91920061 

0.141 

.142 

.143 

.144 

.145 

1.0887071 

1.0895160 

1.0903284 

1.0911444 

1.0919640 

0.91852069 

.91783875 

.91715487 

.91646898 

.91578111 

0.146 

.147 

.148 

.149 

.150 

1.0927873 

1.0936143 

1^0944450 

1.0952795 

1.0961178 

0.91509116 

.91439916 

.91370512 

.91300896 

.91231070 

0.151 

.152 

.153 

.154 

.155 

1.0969599 

1.0978058 

1.0986557 

1.0995095 

1.1003674 

0.91161035 

.91090792 

.91020326 

.90949646 

.908.78737 

0.156 

.157 

.158 

.159 

.160 

1.1012292 

1.1020950 

1.1029650 

1.1038391 

1.1047174 

0.90807618 

.90736280 

.90664708 

.90592913 

.90520888 

0.161 

.162 

.163 

.164 

.165 

1.1055999 

1.1064866 

1.1073777 

1.1082731 

1.1091730 

0.90448633 

.90376151 

.90303426 

.90230468 

.90157261 
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TABLE  I  -  GENERAL  DENSITY  TABLE  -  Continued 
(a)  v  =  1.4  -  Concluded 


2 


5cp 

1 

£ 

Z 

0.166 

.167 

.168 

.169 

.170 

1.1100772 

1.1109859 

1.1118991 

1.1128170 

1.1137394 

0.90083825 

.90010143 

.89936218 

.89862035 

.89787611 

0.171 

.172 

.173 

.174 

.175 

1.1146665 

1.1155983 

1.1165348 

1.1174762 

1.1184225 

0.89712932 

.89638000 

.89562815 

.89487364 

.89411649 

0.176 

.177 

.178 

.179 

.180 

1.1193737 

1.1203298 

1.1212910 

1.1222573 

1.1232287 

0.89335670 

.89259431 

.89182915 

.89106126 

.89029064 

0.181 

.182 

.183 

.184 

.185 

1.1242053 

1.1251871 

1.1261742 

1.1271667 

1.1281647 

0.88951724 

.88874108 

.88796209 

.88718022 

.88639540 

0.186 

.187 

.188 

.189 

.190 

1.1291681 

1.1301771 

1.1311917 

1.1322120 

1.1332381 

0.88560773 

.88481708 

.88402346 

.88322682 

.88242709 

0.191 

.192 

.193 

.194 

.195 

1.1342699 

1.1353077 

1.1363514 

1.1374012 

1.1384571 

0.88162438 

.88081848 

.88000948 

.87919724 

.87838180 

0.196 

.197 

.198 

.199 

.200 

1.1395191 

1.1405873 

1.1416619 

1.1427429 

1.1438304 

0.87756318 

.87674131 

.87591607 

.87508748 

.87425548 

0.201 

.202 

.203 

.204 

.205 

1.1449245 

1.1460252 

1.1471327 

1.1482470 

1.1493681 

0.87342004 

.87258116 

.87173873 

.87089276 

.87004329 

0.206 

.207 

.208 

.209 

.210 

1.1504963 

1.1516316 

1.1527741 

1.1539238 

1.1550810 

L_  _ _ 

0.86919011 

.86833324 

.86747265 

.86660835 

.86574015 

0.211 

.212 

.213 

.214 

.215 

1.1562457 

1.1574179 

1.1585978 

1.1597855 

1.1609811 

0.86486808 

.86399217 

.86311229 

.86222840 

.86134046 

0.216 

.217 

.218 

.219 

.220 

1.1621847 

1.1633964 

1.1646164 

1.1658447 

1.1670816 

0,86044843 

.85955226 

.85865183 

.85774718 

.85683812 

5cp 

1 

2 

Z 

0.221 

.222 

.223 

.224 

.225 

1 .1683271 
1.1695813 
1.1708444 
1.1721164 
1.1733975 

0.85592468 

.85500683 

.85408445 

.85315759 

.85222612 

0.226 

.227 

.228 

.229 

.230 

1.1746880 

1.1759880 

1.1772974 

1.1786165 

1.1799455 

0,85128987 

.85034881 

.84940305 

.84845240 

.84749677 

0.231 

.232 

.233 

.234 

.235 

1.1812846 

1.1826339 

1.1839935 

1.1853636 

1.1867445 

0.84653605 

.84557021 

.84459923 

.84362300 

.84264136 

0.236 

.237 

.238 

.239 

.240 

1.1881362 

1.1895390 

1.1909530 

1.1923785 

1.1938156 

0.84165435 

.84066180 

.83966370 

.83865987 

.83765030 

0.241 

.242 

.243 

.244 

.245 

1.1952645 

1.1967255 

1.1981988 

1.1996846 

1.2011831 

0.83663490 

.83561351 

.83458605 

.83355242 

.83251255 

0.246 

.247 

.248 

.249 

.250 

1.2026946 

1.2042193 

1.2057574 

1.2073092 

1.2088750 

0.83146628 

.83041353 

.82935423 

.82828823 

.82721539 

0.251 

.252 

.253 

.254 

.255 

1.2104551 

1.2120497 

1.2136591 

1.2152836 

1.2169235 

0.82613556 

.82504868 

.82395460 

.82285320 

.82174434 

0.256 

.257 

.258 

.259 

.260 

1.2185792 

1.2202510 

1.2219392 

1.2236442 

1.2253663 

0.82062783 

.81950353 

.81837132 

.81723102 

.81608251 

0.261 

.262 

.263 

.264 

.265 

1.2271059 

1.2288634 

1.2306393 

1.2324338 

1.2342475 

0.81492559 

.81376010 

.81258578 

.81140261 

.81021027 

0.266 

.267 

.268 

.269 

.270 

1.2360808 

1.2379342 

1.2398081 

1.2417032 

1.2436198 

0.80900860 

.80779738 

.80657644 

.80534543 

.80410428 

0.271 

.272 

.273 

.274 

.275 

1.2455586 

1.2475200 

1.2495048 

1.2515134 

1.2535466 

0.80285263 

.80159036 

.80031705 

.79903260 

.79773660 

5-cp 

1 

H 

Z 

0.276 

.277 

.278 

.279 

.280 

1.2556050 

1.2576894 

1.2598004 

1.2619389 

1.2641055 

0.79642881 

.79510887 

.79377654 

.79243139 

.79107321 

0.281 

.282 

.283 

.284 

.285 

1 . 2663012 
1.2685268 
1.2707833  - 
1.2730715 
1.2753927 

0.78970153 

.78831602 

.78691623 

.78550184 

.78407223 

0.286 

.287 

.288 

.289 

.290 

1.2777476 

1.2801376 

1.2825638 

1.2850276 

1.2875299 

0.78262718 

.78116602 

.77968831 

.77819340 

.77668099 

0.291 

.292 

.293 

.294 

.295 

1.2900725 

1.2926566 

1.2952840 

1.2979563 

1.3006755 

0.77515023 

.77360066 

.77203146 

.77044196 

.76883127 

0.296 

.297 

.298 

.299 

.300 

1.3034429 

1.3062609 

1.3091314 

1.3120571 

1.3150400 

0.76719893 

.76554385 

.76386526 

.76216195 

.76043314 

0.301 

.302 

.303 

.304 

.305 

1.3180832 

1.3211890 

1.3243612 

1.3276021 

1.3309163 

0.75867745 

.75689398 

.75508102 

.75323774 

.75136205 

0.306 

.307 

.308 

.309 

.310 

1.3343067 

1.3377787 

1.3413356 

1.3449840 

1.3487277 

0.74945288 

.74750779 

.74552558 

.74350327 

.74143951 

0.311 

.312 

.313 

.314 

.315 

1.3525748 

1.3565299 

1.3606032 

1.3648001 

1.3691338 

0.73933065 

.73717505 

.73496814 

.73270804 

.73038880 

0.316 

.317 

.318 

.319 

.320 

1.3736105 

1.3782479 

1.3830531 

1.3880508 

.1.3932489 

0,72800841 

.72555888 

.72303804 

.72043473 

.71774684 

0.321 

.322 

.323 

.324 

.325 

1.3986833 

1.4043617 

1.4103403 

1.4166227 

1.4233067 

0.71495813 

.71206727 

.70904873 

.70590426 

.70258926 

0.326 

.327 

.328 

.329 

.330 

1.4303770 

1.4380371 

1.4461815 

1.4554567 

1.4650468 

0.69911639 

.69539235 

.69147614 

.68706956 

.68257205 
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TABLE  I 


6qp 

1 

r 

Z 

0.001 

.002 

.003 

.004 

.005 

1.0005006 

1.0010026 

1.0015060 

1.0020107 

1.0025168 

0.99949965 

.99899840 

.99849626 

.99799333 

*.99748952 

0.006 

.007 

.008 

.009 

.010 

1.0030242 

1.0035330 

1.0040432 

1.0045548 

1.0050678 

0.99698492 

.99647944 

.99597308 

.99546585 

.99495775 

0.011 

.012 

.013 

.014 

.015 

1.0055822 

1.0060980 

1.0066152 

1.0071339 

1.0076541 

0.99444879 

.99393896 

.99342827 

.99291663 

.99240404 

0.016 

.017 

.018 

..019 

.020 

1.0081756 

1.0086985 

1.0092230 

1.0097490 

1.0102764 

0.99189070 

.99137651 

.99086129 

.99034513 

.98982813 

0.021 

.022 

.023 

.024 

.025 

1.0108053 
1.0113357 
1.0118676  • 
1.0124010 
1.0129359 

0.98931021 

.98879136 

.98827159 

.98775090 

.98722930 

0.026 

.027 

.028 

.029 

.030 

1.0134723 

1.0140104 

1.0145499 

1.0150910 

1.0156336 

0.98670679 

.98618318 

.98565876 

.98513335 

.98460705 

0.031 

.032 

.033 

.034 

.035 

1.0161778 

1.0167236 

1.0172710 

1.0178200 

1.0183707 

.98407975 

.98355148 

.98302222 

.98249199 

.98196069 

0.036 

.037 

.038 

.039 

.040 

1.0189229 

1.0194767 

1.0200322 

1.0205894 

1.0211482 

0.98142853 

.98089539 

.98036121 

.97982597 

.97928978 

0.041 

.042 

.043 

.044 

.045 

1.0217087 

1.0222708 

1.0228346 

1.0234001 

1.0239673 

0.97875255 

.97821438 

.97767518 

.97713494 

.97659369 

0.046 

.047 

.048 

,049 

.050 

1.0245363 

1.0251070 

1.0256794 

1.0262536 

1.0268295 

0.97605131 

.97550792 

.97496352 

.97441802 

.97387151 

0.051 

.052 

.053 

.054 

.055 

1.0274072 
1.0279867 
1.0285679 
1.0291510 
.  1.0297359 

0.97332392 

.97277523 

.97222556 

.97167471 

.97112279 

-  GENERAL  DENSITY  TABLE  -  Continued 
(b)  7  =  4/3 


2 


6cp 

1 

Z 

Z 

0.056 

.057 

.058 

.059 

.060 

1.0303226 

1.0309112 

1.0315016 

1.0320939 

1.0326880 

0.97056980 

.97001565 

.96946044 

.96890409 

.96834668 

0.061 

.062 

.063 

.064 

.065 

1.0332840 

1.0338819 

1.0344817 

1.0350835 

1.0356872 

0.96778814 

.96722846 

.96666766 

.96610563 

,96554249 

0.066 

.067 

.068 

.069 

.070 

1.0362928 

1.0369003 

1.0375099 

1.0381215 

1.0387350 

0.96497824 

.96441288 

.96384622 

.96327838 

.96270945 

0.071 

.072 

.073 

.074 

.075 

±. 0393506 
1.0399681 
1.0405876 
1.0412093 
1.0418331 

0.96213924 

.96156796 

.96099550 

.96042169 

.95984664 

0.076 

.077 

.078 

.079 

.080 

1.0424589 

1.0430867 

1.0437167 

1.0443488 

1.0449830 

0.95927043 

.95869308 

.95811440 

.95753449 

.95695337 

0.081 

.082 

.083 

.084 

.085 

1.0456194 

1.0462580 

1.0468987 

1.0475416 

1.0481868 

0.95637093 

.95578720 

.95520226 

.95461603 

.95402842 

0.086 

.087 

.088 

.089 

.090 

1.0488341 

1.0494837 

1.0501355 

1.0507896 

1.0514460 

0.95343963 

.95284948 

.95225807 

.95166530 

.95107119 

0.091 

.092 

.093 

.094 

.095 

1.0521047 

1.0527657 

1.0534290 

1.0540947 

1.0547628 

0.95047575 

.94987897 

.94928087 

.94868137 

.94808046 

0.096 

.097 

.098 

.099 

.100 

1.0554332 

1.0561061 

1.0567813 

1.0574590 

1.0581391 

0.94747825 

.94687456 

.94626958 

.94566314 

.94505533 

0.101 

.102 

.103 

.104 

.105 

1.0588217 

1.0595068 

1.0601944 

1.0608846 

1.0615773 

0.94444608 

.94383538 

.94322324 

.94260959 

.94199452 

0.106 

.107 

.108 

.109 

.110 

1.0622726 

1.0629704 

1.0636709 

1.0643740 

1.0650797 

0.94137795 

.94075997 

.94014041 

.93951938 

.93889687 

6cp 

1 

Z 

Z 

0.111 

-.112 

.113 

.114 

.115 

1.0657881 

1.0664991 

1.0672128 

1.0679293 

1.0686486 

0.93827281 

.93764730 

.93702025 

.93639158 

.93576130 

0.116 

.117 

.118 

.119 

.120 

1.069370 6 
1.0700954 
1.0708230 
1.0715534 
1.0722867 

0.93512951 

.93449612 

.93386115 

.93322461 

.93258641 

0.121 

.122 

.123 

.124 

.125 

1.0730229 

1.0737620 

1.0745040 

1.0752489 

1.0759968 

0.93194656 

.93130508 

.93066196 

.93001723 

.92937079 

0.126 

.127 

.128 

.129 

.130 

1.0767477 

1.0775017 

1.0782587 

1.0790187 

1.0797818 

0.92872267 

.92807278 

.92742122 

.92676800 

.92611304 

m 

j||HEE  *9 

plliWir*  *9 

1.0805481 

1.0813175 

1.0820901 

1.0828659 

1.0836449 

0.92545626 

.92479776 

.92413746 

.92347538 

.92281152 

0.136 

.137 

.138 

.139 

.140 

1.0844271 

1.0852126 

1.0860015 

1.0867937 

1.0875893 

0.92214590 

.92147843 

.92080904 

.92013783 

.91946473 

0.141 

.142 

.143 

.144 

.145 

1.0883882 

1.0891906 

1.0899965 

1.0908058 

1.0916187 

0.91878982 

.91811295 

.91743414 

.91675347 

.91607079 

0.146 

.147 

.148 

.149 

.150 

1.0924351 

1.0932551 

1 .0940787 
1.0949059 
1.0957368 

0.91538619 

.91469960 

.91401103 

.91332050 

.91262792 

0.151 

.152 

.153 

.154 

.155 

1.0965714 

1.0974098 

1.0982520 

1.0990979 

1.0999477 

0.91193332 

.91123662 

.91053784 

.90983706 

.90913413 

0.156 

.157 

.158 

.159 

.160 

1.1008014 

1.1016591 

1.1025207 

1.1033863 

1.1042559 

0.90842908 

.90772182 

.90701245 

.90630090 

.90558719 

0.161 

.162 

.163 

.164 

.165 

1.1051296 

1.1060074 

1.1068894 

1.1077756 

1.1086660 

0.90487125 

.90415308 

.90343263 

.90270990 

.90198491 
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TABLE  I  -  GENERAL  DENSITY  TABLE  -  Concluded 


(b)  r  =  4/3  -  Concluded 
2 


6cp 

1 

Z 

Z 

0.166 

,167 

.168 

.169 

.170 

1.1095607 

1.1104597 

1.1113630 

1.1122708 

1.1131830 

0.90125759 

.90052795 

.89979602 

.89906163 

.89832489 

0.171 

.172 

.173 

.174 

.175 

1.1140997 

1.1150210 

1 . 1159468 
1.1168773 
1.1178125 

0.89758574 

.89684410 

.89610006 

,89535350 
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Figure  5.  -  Mean  stream  surfaces  for  inlet  stage  of  axial-flcw  compressor 
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(a)  surface  with  qp  and  z  as 
independent  variables. 


(b)  S1  surface  with  r  and  cp  as 
independent  variables. 


(c)  S2  surface  with  r  and  z  as 
independent  variables. 
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Figure  7.  -  Elements  of  stream  sheet. 


12E 


89 


MCA  TF  2604 


Figure  8.  -  Orthogonal  coordinates  for  surface  of  revolution. 
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Figure  10.  -  Grid  system  and  "boundary  conditions  for  general  surface 

(elliptic  case). 
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Figure  11.  -  Grid  system  and  boundary  conditions  for  general  S2  surface  (elliptic  case 


